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ABSTRACT 


This thesis presents the use of a general purpose non- 
linear optimization program in the preliminary stage of ship 
design for the selection of a propeller based on methodical 
series propeller test data. The propeller series utilized is 
the well-known Wageningen B-Series. Three (3) "Design Cases", 
representing the thrust, power and matching approaches to 
powering problems, are formulated as FORTRAN subprogram 
analysis codes for solution by the synthesis/optimization 
program COPES/CONMIN. Designer constraints considered are: 

1) diameter limitation : 


2) cavitation limit on expanded area ratio using 
Keller's criterion 


3) strength requirement determined by an empirical 
relation and by a method developed by Schoenherr 
with modifications by the author. 
Objective functions considered are maximized open water 
efficiency and minimized propeller blade weight. Optimized 


solutions to specific problems previously presented by other 


authors are obtained and results are compared. 
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A. BACKGROUND 

The ship design process, in its most rudimentary form, 
has been formulated and tracked by the utilization of the 
classical design spiral (see Figure 1.1). The design follows 
a convergent helical path past each major milestone "spoke" 
until, after numerous iterative cycles, the final configura- 
tion is "centered" upon. Whether one attempts to segregate 
the principal phases of Preliminary, Advanced and Contract 
Design into separate spirals or combine these phases in 
series along the entire path to the center, it is not long 
before the designer's roughed-out sketches give way to serious 
Paumoer Crunching", specifically that of propulsion power 
estimation. 

To estimate the power required to drive the ship through 
the water at its design speed, a decision must first be made 
as to what type of propulsor (1.e., propeller, water jet, 
paddle wheel, etc.) will be used. For the average case, and 
for the discussion that follows, the marine propeller is 
chosen to be the propulsion device. Since 

a ship propeller may be regarded as a transducer 

that converts the rotational power transmitted 

through the shaft into the translational power to 

juoecl! the ship, [Ref. 18: p. 10) 
the selection and design of this device is obviously an 
important factor in the eventual size (weight and power) of 


the ship's propulsion plant. While hydrodynamicists provide 
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a myriad of theories and techniques to generate a "Custom 
built" (i.e., wake adapted) propeller for the ship under 
consideration, their expertise is usually not required in 
the early stages of preliminary design simply because the 
design has not been refined enough beyond gross estimates. 
At this stage, the designer strives to formulate what is 
possible based on previous experience. For preliminary 
power estimation, previous propeller designs (i.e., "stock" 
propellers) and results from methodical series of model 
propellers are analyzed by the designer in order to select 
the "best" available propeller under various conditions posed 
by the problem under consideration. Three examples of typi- 
cal problems encountered in preliminary ship design are: 

1) Given the ship's effective horsepower at a specific 
speed and estimates of hull performance parameters, which 
propeller, as determined by certain principal characteristics, 
will require the least amount of delivered power from the 
Propulsion plant? 

2) Given the delivered power from a specific propulsion 
system in terms of torque and revolution rate at the pro- 
peller/shaft interface and estimates of hull performance 
parameters, which propeller will generate the largest effec- 
tive horsepower and speed parameters? 

3) Given a ship's effective horsepower and speed, various 
hull performance parameters, and the propulsion plant's 


delivered power characteristics, which propeller will "match" 


iG 





these requirements at a minimum amount of weight for a speci- 
fied material? 

(Author's Note: For the sake of brevity, the three 
selection problems just cited will, henceforth, be referred 
to as "Design Case No. 1", "Design Case No. 2" and "Design 
Case No. 3", respectively.) 

For this study, the methodical propeller series method 
is viewed as the designer's choice for preliminary powering 
analysis. One of the most-widely used methodical series data 
on model propellers is the Wageningen B-Screw Series. 
Initially, the results of the series were presented as tabu- 
lations of non-dimensional thrust and torque coefficients 


(Ky and K. respectively) versus the non-dimensional advance 


Q 
ratio (J) for analytical work anc as the familiar "Bp-6" and 
"Bu-é" diagrams for design purposes. As "trial & error" 
design methods performed by hand in all engineering disci- 
plines gradually transcended to numerical manipulation by 
the modern digital computer, the necessity for the adaption 
of the Series results to a format suitable for use in computer- 
aided design methods became obvious. This was accomplished 
through multiple regression analysis of the original open- 
Water test data of the 120 propeller models in the Series 
and presented in the form of polynomial expressions for "Kye 
and "Ko" (REES= -L,21). 

The adaptation of the Wageningen B-Screw Series polynomials 


to various types of propeller selection problems formulated 


ale 





for computer solution has been implemented recently by two 
Piene aoe torancadtyllOu |Reft. 3) and, Of late,, Markussen 
[Ref. 4] presented different propeller selection problems and 
proposed different schemes for computer-aided "optimized" 
solutions. In short, specific expressions for the con- 
straints imposed and the objective (optimality condition) 

to be maximized, expressed in terms of a number of design 
variables and parameters, were developed. Then, each system 
of equations was solved by a Newton-Raphson method to give 

a solution set of the design variables which maximized the 
objective and met all constraints. 

Rather than formulating and coding a different optimiza- 
tion scheme each time a propeller selection problem presents 
a different combination and number of design parameters, 
variables and constraints, a better approach would involve 
formulating the problem (constraints and objective function) 
once in terms of all design parameters and variables and 
utilizing a general purpose optimization scheme which can 
handle any combination and number of constraints and design 
variables. This alternative certainly allows the designer 
more flexibility in solving his problem. Moreover, it elimin- 
ates repetitive coding and debugging associated with the 
implementation of a computer-sided solution for each particu- 


lar design problem. 
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B. PROBLEM STATEMENT 

The problem, then, is that the previously cited computer- 
aided "optimized" solutions to the propeller selection problem 
are not broad enough in capability to handle variations in 
the problem formulation. The objective of this thesis is to 
apply an available general purpose optimization computer 
code to the solution of various propeller selection problems 
encountered in Preliminary Ship Design in order to enhance 


the flexibility of the selection procedure. 


fe. SCOPE 

To achieve the stated objective, the general purpose 
non-linear optimization code CONMIN [Refs. 5,6] together with 
the engineering synthesis code COPES [Ref. 7] (hereafter 
referred to collectively as COPES/CONMIN) is utilized in the 
solution of the three previously cited preliminary design 
propeller selection problems. Using the Wageningen B-Screw 
Series propeller characteristics expressed in polynomial 
expressions of various design variables, three "analysis" 
codes, required by COPES/CONMIN, are developed in such a way 
that various combinations of design variables and constraints 
are used, thereby demonstrating the applicability of the 
COPES/CONMIN optimization program in the solution of propeller 


selection problems. 


D. THESIS ORGANIZATION 
The remainder of the thesis is organized in the following 


manne:xr . 


i 





Chapter II presents a short description of the optimiza- 
tion problem in general terms and a follow-on discussion of 
the COPES/CONMIN optimization program and the Potomac cal 
techniques employed therein. 

Chapter III introduces definitions and concepts applica- 
ble to the propeller selection problem. A subsequent dis- 
cussion on the Wageningen B-Screw Series is followed by final 
comments on constraints imposed on the propeller selection 
preplem. 

Chapter IV presents the formulation of the propeller 
selection problem as a design optimization problem which 
can be solved using COPES/CONMIN. 

Chapter V discusses the background, formulation and 
programming utilized in estimating a propeller blade's weight 
for subsequent consideration as an objective function. 

Chapter VI reviews the author's modifications to the 
propeller strength analysis developed by Schoenherr [Ref. 8] 
in the early 1960's for the American Bureau of Shipping. A 
subsequent discussion on the programming details of FORTRAN 
codes, which are utilized for the determination of adequate 
propeller blade strength, completes the chapter. 

Chapter VII reviews the formulation and programming for 
the analysis code which is used in solving propeller selection 
problems represented by Design Case No. 1. Sample solutions 
are presented and compared to those presented previously by 


other authors. 
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Chapters VIII and IX consider Design Case No. 2 and 
Design Case No. 3 selection problems, respectively, in a 
Similar fashion to Chapter VII. 

Chapter X, the final chapter, presents the author's 
conclusions and recommendations. 

As a final note, all computer coding presented in this 
thesis is done in FORTRAN IV, the language used by COPES/ 
CONMIN. For the reader's convenience, Appendix A provides 
a cross-reference of the symbols presented throughout the 
thesis to appropriate FORTRAN variable names appearing in 


the author's codes. 
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Pee LN TRODUCTION 

The purpose of this chapter 1s to introduce definitions 
and concepts used in the formulation and solution of the 
general optimization problem. Then, a short discussion on 
the theory and implementation details of COPES/CONMIN is 
presented. 

For further study on the theory and methods of optimiza- 
tion, the reader is directed to the texts by Fox [Ref. 9], 


Fiacco and McCormick [Ref. 10], and Himmelblau [Ref. 11]. 


B. DEFINITIONS 

Before discussing the techniques of optimization and 
their application to engineering problems, some preliminary 
definitions of basic terminology should be stated. Terms 
which have relevant significance are: 

1) Parameters--The numerical quantities for which values 
are assigned to produce a design are called parameters. From 
this, it follows that a design may be specified by a vector 
D containing "p" components, each of which is associated with 


a parameter. That is: 


a é p22uL} 
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However, in a design process, the parameters are determined 
by some logical procedure through analysis of some kind. 
Some might take on fixed values to become "“preassigned" 
parameters. Interrelationships among other parameters might 
exist so that only some of the parameters are changed when 
one design is compared to another. This consequence leads 
to the definition of "design variable". 

2) Design Variables--The parameters for which values are 
chosen in some fashion to produce a de3icn are called design 
variables. They represent an ordered collection of components 
which is a subset of the design vector D. This subset is 
unique in that its components are "variable", i.e., they may 
take on different values in the design process. Having 
"preassigned" or fixed some of the design's parameters and 
only allowing the remaining "design variables" to change, leads 
to the conclusion that a design is now uniquely specified by 
a vector X containing "n" components (n < p), each of which 


is associated with a design variable. That is: 


\>] 
X = os ry 


3) Objective Function--The computable function of all or 
some of the deSign's preassigned parameters and/or design 
variables, with respect to which the design is to be opti- 


mized, is called the objective function. Single valued in 


Ze 





Guantitative terms, the objective function's minimum or maxi- 
mum value represents the "best" obtainable or "optimized" 
design. It is expressed as F(D) to show its dependence on 
the design's parameters. But, since a design can be uniquely 
defined by X alone, then clearly F(X) suffices as an expression 
For the objective function. 

4) Constraints--Restrictions on the design which must be 
satisfied in order to produce an acceptable design are 
called constraints. A constraint may be classified as a 
"Side" or a "behavior" constraint. A side constraint re- 
stricts or bounds the range of the design for reasons other 
than direct consideration of performance. The side constraint 


on the "i"th design variable may be expressed as: 


lower < X. eX 


Xx 
aE = ii a 


A constraint derived from those performance or behavior 
requirements that are explicitly considered is called a be- 
havior constraint. Most often, it appears as a computable 
functional relation involving the design's parameters, both 
preassigned and variable alike. The relation may be an in- 
equality so that the "j"th of "m" inequality constraints 


can be expressed as: 


G, (D) <——e ee lor a4) 
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Alternatively, the relation may be an equality on the "k"th 


of "2" equality constraints expressed as: 


Hy (D) = 0 ee Oy ouetony (275) 


Of noteworthy importance here is that, as before, if 
some of the design's parameters are preassSigned, then the 
resulting design is that defined by X which contains only 
the parameters that can be varied in the design process, 
1.e., the design variables. Therefore, constraints imposed 


upon the design may be expressed under one equation as: 


ne 7, 0 Gil. 7m 
Hy, (X) = 0 KS a heat (27.16) 
ae 2 on = axe Se iets: abete-c. 2 oa 
a — ir 1 


Pecinal form for constraints is that of the discrete-valued 
design variable. 

5) Feasible Design--A design in which specified constraints 
are satisfied is called a feasible or "acceptable" design. 

6)  Infeasible Design--A deSign in which constraints are 


violated is called an infeasible or "unacceptable" design. 


C. PROBLEM STATEMENT 
If one presupposes that a range of designs exists within 


a selected design concept, then it follows that different 
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methodologies also exist by which one may choose the param- 
eters which describe the design. One such method is optimi- 
zation where parameters are chosen in a way that the sehen 
will satisfy all of the limitations and restrictions imposed 
upon it and will be "best" in some sense. In view of the 
foregoing definitions, optimization is then a selection method 
applied to a design problem by which an objective function 
F(D) is minimized to produce an acceptable design which satis- 
fies a certain set of requirements called constraints. 
Formulated mathematically, the general, non-linear, 


constrained optimization problem may be stated under one 


equation as: 


Minimize: F(D) = OBJ 
Subject to: G, (D) < Bae. leer ail (25 7) 
Hy, (D) ono eee lames: 
ylower <0) Xone eee at yeah awa er 
iE —- 1-1 


Again, aS pointed out in the previous section, the design 
May be uniquely defined by just its design variables as 
specified by X when some parameters are preassigned. Thus, 
the general, non-linear, constrained optimization problem 


can now be stated under one equation as: 
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Minimize: F(X) = OBJ 


Subgects to. 8 G(X) <9 0 pee ly tpn (2e8)) 
a eC koge Cl corneas 
ee Seteeweey et = 1,.5.,;n 


Solutions methods for this optimization problem are 
abundant. Those pertaining to the linear and quadratic 
optimization problems involving a few deSign variables are 
most often presented in graphical or analytic form, although 
numerical schemes are, by no means, a dormant form. Struc- 
tural and thermal problem solutions are most prevalent. How- 
ever, as the optimization problem becomes more complex in 
terms of non-linear relationships among an increasing number 
of design variables and of an increased number of design 
constraints, numerical or mathematical programming techniques 
dominate the solution methods. 

To limit the scope of this discussion, only the numerical 
technigues relevant to COPES/CONMIN will be considered. For 
more background on optimization techniques and applications, 
the reader is directed to a recent paper by Vanderplaats 
(Ref. 12] which presents a concise, but thorough, qualitative 
review of optimization. Although this paper deals exclusively 
with the application of design optimization to structural 
problems, it also contains a very extensive and current list 
of references on general techniques and applications of 


optimization. 
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De COrESs/ CONMIN 

As previously stated in Chapter I, COPES/CONMIN is the 
collective acronym for the FORTRAN program utilizing the opti- 
mization code CONMIN and the synthesis code COPES. COPES 
stands for COntrol Program for Engineering Synthesis; CONMIN 
is an acronym for CONstrained function MINimization. 

1. CONMIN 

CONMIN is a FORTRAN program, in subroutine form, which 

solves the general non-linear constrained optimization prob- 


lem as stated: 


Minimize: F(X) = OBJ 
Subject to: G(X) a 0 een (229) 
OWE. ee ec s | e 
a = = L 


Equation (2.9) applies to the entire statement. Observe 
that equation (2.9) differs from equation (2.8) in that the 
equality constraint set, given by HX) = 0 ersenotmspee1 = led . 
This is because the version of COPES/CONMIN used in this 
study does not consider these types of constraints. However, 
this will not pose any difficulty in solving the propeller 
selection problems previously cited. 

Again, F(X) is the objective function (OBJ). The 
vector X contains the "n" design variables (NDV). G(X) 


are the "m" inequality behavior constraints (NCON) imposed 


ylower upper 
1 


on the optimization problem; and X are the 


Zo 





respective lower and upper side constraints which bound the 
"design space" over which F(X) and G, (x) are defined. As 
functional relationships involving xX, F(X) and G(X) may be 
mmplicit Or Explicit, but, in any event, must be continuous 
and have finite numerical values. 

When the inequality condition of equation (2.9) is 
Metesatictied, i.e., G, (X) > 0 for any constraint, the con- 
straint is said to be violated. If the equality condition 
Mmeormet, 1.€., G, (x) =Uetomeamywconstraint, the constraint 
is said to be active. And, finally, if the inequality condi- 
mioneis satisfied, i.e., G; (x) 0) for any Constraine, that 
constraint is termed inactive. Any design, defined by Xx, 
which satisfies the inequalities of equation (2.9) is desig- 
nated as a feasible design. Likewise, any one which vio- 
lates these inequalities is termed an infeasible design. 

The feasible design with the minimum objective function value, 
often referred to as the "minimum feasible design", will, 
therefore, be the optimum design. 

During the optimization process, CONMIN employs the 
Fletcher-Reeves algorithm [Ref. 13] for locally unconstrained 
problems, and Zoutendijk's method of deasible directions [Refs. 
14,15] for locally constrained problems, in a numerical pro- 
cedure which attempts to minimize the objective function, 
F(X) = OBJ, until one or more of the constraints, G(X), 
becomes active. The numerical search procedure begins with 


an initial X vector which may or may not specify a feasible 
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design. Modifications are included in CONMIN so that, if 
the initial design is infeasible, a feasible solution will 
be obtained with minimal increase in F(X). By iteratively 
updating the design vector X by the following relation: 
ex oS () Ganon 

the optimization process continues by following the con- 
Seraint boundaries in a direction of search § so that the 
value of F(X) decreases with each iteration q. The scalar 
a* defines the distance of travel in the direction of search 
S. The process terminates when a vector X is found such 
that no further decrease in F(X) can be made. The vector X 
1s considered to be optimal and, at least, a local minimum. 

CONMIN can be used alone as a subroutine in any 
FORTRAN program where numerical optimization is desired. 
However, in order to make the optimization process more "user- 
friendly", CONMIN has been coupled to COPES in order to 
Simplify its application to various types of problems. 
Further information on CONMIN can be found in previously 
cited references [5] and [6]. 

Zoe COLES 

COPES is a FORTRAN program that provides automated 
design and trade-off capability to the design engineer. It 
utilizes the optimizer CONMIN to provide the following six 


specific capabilities: 
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1) simple analysis 

2) Opeimizatilon 

3) sensitivity analysis 

4) two variable function space analysis 

5) optimum sensitivity 

6) optimization using approximation techniques 
During the execution of COPES, say for optimization, three 
principal tasks are performed: 

1) data management on the design variables and constraints 
through location assignments in a FORTRAN common block called 
GLOBCM. 

2) decision process control on the attainment of an optimal 
design vector X through multiple calls to the optimizer 
until a minimum or maximum value of OBJ is achieved and all 
G, (X) are satisfied. 

3) evaluation of OBJ and G., (X) at each X@ and atl when 
ICALC = 2 through multiple calls to the user-provided analysis 
subprogram, SUBROUTINE ANALIZ. 

For the application under consideration in this study, only 
the optimization capability will be used. Therefore, further 
elaboration on the other capabilities is not warranted. 
Reference [7] is the user's manual for COPES/CONMIN. 
Details on the mechanics of user implementation are presented 
with subsequent illustration by example. The reader is, 
therefore, encouraged to familiarize himself with the refer- 
ence. However, at this point, it is sufficient to be aware 


of the fact that a user of COPES/CONMIN is required to: 
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1) provide a FORTRAN subroutine called ANALIZ which per- 
forms the input of preassigned parameters, the evaluation of 
the objective function and constraints during the analysis 
phase of the optimization search and the output of the 
results. 

2) provide an assembled deck of control cards required 


iy COPES. 


Be CONCLUDING NOTE 

The field of optimization is both extensive and complex 
and, therefore, the foregoing presentation is, by no means, 
complete in every detail. However, it is felt that the pre- 
Be ding overview, in conjunction with the cited references, 
covers the necessary prerequisites that will enable the reader 
to follow the application of COPES/CONMIN to the various 


propeller selection problems in the chapters that follow. 
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III. POWERING, PERFORMANCE AND PROPELLERS 


A. INTRODUCTION 

The purpose of this chapter is to present an overview 
of the terminology and concepts that pertain to ship propul- 
sion, propeller selection and the use of model propeller 
test data. Initially, fundamental definitions used in ship 
powering problems are presented. This is followed by a 
discussion of the "classic" types of propeller design/ 
selection problems encountered by the naval architect and 
Marine and naval engineers. Propeller model testing and 
propeller performance characteristics are reviewed next. 
The chapter is completed with a discussion of the Wageningen 
B-Screw Series. 

The goal here is brevity. The reader is, therefore, 
encouraged to investigate the references cited for further 


details. 


mee DEFINITIONS 
Some fundamental terms associated with most propeller 
design/selection problems are: 


1) Effective Horsepower (P_)--power required to tow the 


13; 


"bare" hull (without propeller; rudder and appendage allowance 





assumed included) that generates a given resistance (Ry) at 
a given speed (V). It is determined by: 
Re 
ee 
PE ~ 550 Copel 


SS 





2) Thrust Horsepower (P,,) -~power delivered to water by 
a propeller developing a thrust force (T) and moving at a 
speed of advance (V.) without the influence of a hull form 


ahead of it. Pa is determined by: 





eso (3.2) 


3) Delivered Horsepower (P,) --power delivered by shaft 
to propeller, normally specified at the outboard side of the 
stern tube. Qq 1s the torque delivered to the propeller; 
a 1s the revolution rate of the shaft and, consequently, 


the propeller. P, is determined by: 


D 
Zi (Ome 
- s P 
Py a samen oi) 
4) Shaft Horsepower (P..) --power delivered to the inboard 


Side of the stern tube having a transmission efficiency of 


Ne: P. 1s determined by 


Po = —— (334) 


5) Brake Horsepower (P.,) --power delivered by the prime 
mover at connection flange to the power train. While P, is 
normally associated with the prime mover's rated power at 


this connection (BHP), it can also be specified from the 


power train/propeller side as: 
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Ree 2s. (3.5) 
: Dena 
where Np and n, are, respectively, the bearing system and 
reduction gear transmission efficiencies. 
6) Thrust deduction factor (l-td)--ratio of the tow 
resistance (Ro) to the thrust (T) provided by the propeller. 


Tt is determined by 
(1-td) = —- (26) 


7) Wake Factor (Taylor's) (l-wt)--ratio of the speed of 


advance (Vq) to the ship's speed (V). It is given by: 
V 
Sh ts 
(l=we) = 7 (src) 
8) Advance Ratio (J)--a non-dimensional value, associated 


with propeller test data presentation (see figure (3.1), 


given by the following relation: 


V (l-wt) A 
eS (3.6) 
ny D Ny D 
9) mirust CocftfLlcient (Kn) --a non-dimensional value asso- 


Clated with the thrust force (T) developed by a propeller 


Of diameter D which is turning at a rate n_ and operating 


P 
in a fluid of density p. It is defined by the following 


expression: 


3D 





10) Torque Coefficient (K,)--a non-dimensional value 
associated with the torque (Q,,) absorbed by a propeller of 


diameter D which is turning at a rate n. and operating in 


e 
a fluid of density op. It is defined by the following 


expression: 


(3210) 


ll) Open Water Efficiency (n,) --the rat LO von Pr to Py 
for a propeller in open water conditions, with a uniform 
inflow velocity field at a speed of advance Vas feels ne x— 


pressed as: 





No = = ~ a 7 - a (3.11) 
D p “Ss ~~ Q 


12) Hull Efficiency (n,,)--a ratio of work done on the 


H 
ship to that done by the propeller expressed as: 


— Peo | Rp VY (ata) 
irre) Ty. (L=wt) Sold) 





13) Relative Rotative Efficiency (np) --the ratio of the 
actual, behind-hull efficiency to the open water efficiency. 


The value of Np does not, in general, depart from the value 
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Sel,0., §—Most Often, Ne varies between 0.95 and 1.0 for twin- 


Peiee wes ps aie oetweenm 70 and 1.1 for single screw ships. 


Applicable units for the terms in the expressions above 


are: 

1) horsepower (hp) --P., Pre Por Ph and Pp 

2) pounds (lbf)--R,, ie 

3) feet/second (ft/sec)--V, Va 

4) foot-pounds (ft-lbf)--Q., Qp 

5) feet (ft)--D 

6) revolutions/second (rps) --n, 

7) yrevolutions/minute (rpm) --N,, = np/60.0 

Bhe quantities T, V., D, Qp and nN, are obtained from the 

propeller test data results. The quantities Ro and V are 


specified from the design point on the R-V curve for the hull 
under study. The quantity op iS a property of the fluid in 
which the hull and propeller operate. And, finally, Npr Ne 
and Ng are characteristics of the bearing, gear and stern 
tube systems. In preliminary design studies, nominal values, 
based on previous designs, are usually assumed unless, of 


course, these systems have been selected and actual values 


can be specified. 


C. POWERING CONCEPTS 
1. Basic Relations 
Simply stated, the fundamental powering relationship 


to be solved in ship propulsion and powering problems is: 


oe 





(hea 


pcan eee) > p (3.13) 


Utilizing the definitions just presented, equation (3.13) 


can be rewritten as: 


alee posta), PAL o: cee? (3.14) 
550 (i-wt) ° "R No 550 
Rearranging terms of equation (3.14) gives: 
ewes ee see ote 
(1-td) 550 NR No 550 


And, finally, when substitutions are made, equation (3.15) 


becomes: 
A ~ TOU-wt)vV _ em Qs Np a 
550 550 = iR G 550 


Mmemetions (3.14), (3.15), and (3.16) provide the basis for 
different approaches to the solution of a typical powering 
problem. More background and information on the definitions 
and equations presented above may be found in Chapter VI, 
Sections 10-16 in the text by Comstock [Ref. 16] and O'Brien's 
book [Ref. 17]. 
g. Approaches to the Powering Problem 

From equation (3.16), three types of propeller selec- 
tion problems are discernible. In the first instance, the 
propeller thrust T and the propeller's speed of advance V 


A 
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are taken as known quantities. The fact that T is known 
substantiates the "Thrust Approach" nomenclature given to 
this type of selection problem. In the preliminary (or, in 
some circles, conceptual) ship design phase, the specifica- 
tion of T is based upon the requirement imposed by the 


resistance of the ship (R at its design speed (V) (or, the 


é 
effective horsepower (P.) at V) and estimates of wt and td 

in the absence of wake surveys and self-propulsion data from 
model tests. Essentially, the thrust delivered by a selected 
propeller must provide, at least, the thrust required for the 
ship hull under study. The objective in the "Thrust Approach" 
selection problem is to determine, by logical means, the 


appropriate values of Q. and n, when the open water efficiency 


e 
(n,) 1s set by the selected propeller and its performance 
characteristics. 

In the second instance, the delivered torque (Q.) and 
the propeller shaft speed (n,) are taken to be known. The 
"Power Approach" nomenclature iS given to propeller selection 


problems of this type because P. is known. Here, with the 


D 
Shaft and propeller speeds being equal, the torque absorbed 

by the propeller (O,) must be, at least, equivalent to the 
delivered torque (Q,) . The corresponding objective in the 
"Power Approach" selection problem is to determine, by logical 
means, the expected ship speed (V) (or, the speed of advance 


(V and the associated thrust (T) that can be developed 


a)? 
when the open water efficiency Sih 1S, again, set by the 


selected propeller and its performance characteristics. 


oo 





The final, and most familiar, types of propeller 


selection problem occurs when T, V, or V, Qc and n 


a are all 


2 
known. From equation (3.16), the open water efficiency cn) 
is now established as a requirement to be met. The objective 
is, simply, to select a propeller whose open water efficiency 
(no). developed thrust and absorbed torque are equivalent to 
or "match" the requirements imposed. Obviously, this approach 
on the selection problem has been designated as a "matching 
Eeoblem” . 


The reader is directed to the paper by Vassilopoulos 


[fRef. 18] for further information. 


D. PROPELLER PERFORMANCE CHARACTERISTICS 

Up until the late 1950's, much of the knowledge about 
the performance of propellers has been gained from experience 
with models. To study the relationships governing their 
behavior, a model propeller is built and run in a towing 
mame witnout any hull ahead of it. This is done by running 
the propeller on a long shaft projecting well ahead of a 
narrow, hydrodynamically shaped pod or "propeller boat" which 
contains the driving mechanism and recording apparatus and 
is attached to the towing carriage. The propeller advances 
into undisturbed fluid (usually water of density p and kine- 
matic viscosity v) so that the speed of advance (V.) is known 
and the flow into the "disc" swept by the turning blades is 
uniform. For the model propeller of diameter (D) under test, 


meadings of thrust (T), torque (Qo) and shaft revolutions (n,) 
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) 


are recorded over a range of values for speed of advance (V 


in this “open water" condition. : 

Using the laws of similitude, the collected data is 
reduced and scaled appropriately into the familiar functional 
relationships between the advance ratio (J) and the non- 
dimensional coefficients of propeller performance. These 
coefficients or performance characteristics, defined pre- 
viously, are: 

1) Thrust Coefficient (Ky) 

Zee loeaque Coefficient (Ko) 

3) Open Water Efficiency (n,) 
Figure (3.1) graphically depicts the relationship between J 
and Ka and Ka derived from test data for a propeller defined 
by a specific expanded area ratio (AD/A,) s pitch-diameter 
ratio (P/D), number of biades (Z) and thickness-to-chord 
feat do (t/c). 

Definitions of these terms with graphical illustrations 
pertaining to various aspects of propeller geometry can be 
found in Section 15 of references [16], [17] and in van Manen's 
SiblLication [{[Ref. 19]. 

More recently, highly analytical theories (lifting line, 
Meaitied lifting line, lifting surface, etc.) for use with 
high-speed digital computers have been formulated and subse- 
quently used in "modeling", in a mathematical sense, the 


propeller and its behavior in the "wake adapted" (or, behind 


hull) condition as well as the "open water" condition. 
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Additional benefits derived from this approach to propeller 
performance analysis include: 

1) determination of blade section profiles along the 
propeller's blade radius (R) to achieve uniform lift and 
internal stress distributions; 

2) computation of "off-design" performance characteris- 
ees in all quadrants; 

3) subsequent determination of hull surface forces, bearing 
loads and spindle torques induced by the propeller; 

4) prediction of steady and unsteady stress distributions 
in the propeller blade using the finite element method on 
the blade of the propeller under study. 
Obviously, this approach to propeller performance analysis 
Serves to; 

1) eliminate the time-consuming and expensive model 
construction and testing of propellers in tow tanks and 
Cavitation tunnels; 

2) eliminate the "scaling" discrepancies which inhibit 
the reliability of design charts and model propeller data; 

3) eliminate those design charts altogether. 

As in the case with model experiments, however, the ultimate 
objective remains the same, i.e., establishing the performance 
T Ko and as 


functions of J. Having these relationships enables the ship 


characteristics of the propeller in terms of K 


designer to proceed in solving the power equation (equation 


(3.13)) through any of the approaches previously discussed. 
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E. THE WAGENINGEN B-SCREW SERIES 
l. Background 
The model test data of the Wageningen B-Screw Series 
have been selected for use in the powering problems to be 
solved utilizing COPES/CONMIN. The choice was driven by the 
following considerations: 

1) the Series is widely known and, despite its growing 
obsolescence, is still used in preliminary ship design 
studies. 

2) the availability of previous investigations [Refs. 3,4] 
which utilized the series, for comparative analysis of optimi- 
zation results. 

3) the applicability of the polynomial expressions for K 


as 


and K, to computer-aided analysis. 


Q 

The Series tests were conducted from 1940 through 1960 and, 

therefore, represent propeller designs (principally naval 

and merchant applications) and design philosophy of that era. 
Specifically, the Series consists of 120 model 

propellers. As is customary in methodical or systematic 

model propeller series testing, the number of blades (2), 

expanded area ratio (AL/A) and pitch-diameter ratio (P/D) are 

varied systematically, while the blade outline, the profile 

of the blade's cross section along the blade radius, blade 

cross section maximum thickness (t), blade section chord 


length (c), diameter (D) and propeller hub-to-diameter ratio 


(d/D) were kept constant for given values of Ay/Ay and Z. 


=4 





Table (I) summarizes the variations in Z and AD/A, for each 
set of model propellers having pitch-diameter ratios (P/D) 


Sete, 0.0, O.8, 1.0, 1.2 and 1.4. 


TABLE 1 


Summary of the Wageningen B-Screw Series 


Blade number Z fade area ratio 





2. Series Results 
The test results of the Wageningen B-Screw Series were 


Originally presented in the form of Bp-s, Bu-6 and K K 


nee QO’ 
and -J diagrams. As stated in Chapter I, multiple regression 


analysis was performaed (again, [Refs. 1,2]) on the results 


and Ky: The 


open water efficiency Cie acuaEUneceton Of J, follows £rom 


to produce the polynomial expressions for K 


equation (3.11). The correction for "scale effects" was 
achieved by using Lerb's method of equivalent profiles 
[Ref. 20]. Although Triantafyllou's thesis [Ref. 21] sug- 
gests an improved method for scale correction, the results 


of Reference [2] will be used in this study. 
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In propeller selection problems which use this 


Series, values for Kn and Ko are defined as: 


an 
| 


£1 5,P/D,Ap/AG 1 2,t*/CS ap) 


L 
(32) 
= * 
NG £5(J,P/D,Ap/A5 12, /C 7p) 
To compute K,, and Ko: the following equations are used: 
Kn = Kn 5 AKn 
ale dese 
= ' + 
NG Ko AKo 


The polynomial expressions found in Tables (5) and (6) of 
Reference [2] are then used to evaluate to components K,, 


K! 


Q! AK a and AK... 


Q 
Table (5) in Reference [2] lists the coefficients used 


in the polynomial expressions for K/, and Ki} at an equivalent 


a Q 
Reynolds number (Rn 75R) of 2 iO 2 It is defined as: 
/ 2 2 
Me TSR a Ep a 
ea SIR v : 
where: 
v = kinematic viscosity of the fluid (£t7/sec) ; 
C scp = blade section chord length at 3/4 propeller 


radius (./75R) in feet (ft). 
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fmaccoumte for other effects", coefficients AK and AKo are 
introduced. Table (6) in Reference [2] lists the coeffi- 
cients used in the polynomial expressions for these coeffi- 
cients. "Other effects" include the operation of the 
propeller at an equivalent Reynolds Number different from 

2 aa Also, variations in other parameters which define 
the propeller's geometry, specifically t/c values different 
from the ones fixed by the Wageningen propellers, are taken 
into account by corrections to the equivalent Reynolds number. 
By keeping the blade section's chord length (c) at the value 
of the Wageningen propeller, a change in a blade section's 
maximum thickness from the standard one defined by the Series 
(t ) to one preferred in the selection (t* ) produces a 


new equivalent Reynolds number (Rn* ) given by: 


oR 


P 1+2(t/fe 75p) 
RN scp = exp [4.6082 “ i Ita (t*7o a (Qn Rn 35p74-6052) 


(3.2 0) 
where: 
Boe = the new equivalent Reynolds number; 
hen Cc) 75R = new equivalent t/c at 3/4 propeller 
: radius; 
Rn = the Reynolds number computed by equation 
.75R 
ere Ore 
(ec) 75R = standard equivalent t/fc at 3/4 propeller 


radius for Wageningen propellers. 
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For the Wageningen B-Screw Series, the standard equivalent 


lye is given by: 


MOnOnoo. = OMOOL252Z) Z 


t/e az5R = 2.073 A,/A, (3.21) 


Further details on blade section geometry will be addressed 
in Chapter V. Reference [2] contains background and other 
information on the equations above. 
fe) Lamitations on Series Data 
In utilizing the Wageningen B-Screw Series in any 
propeller selection problem, the following restrictions 
apply to the Series data: 
1) Number of Propeller Blades (Z)--The Series considers 
only propellers with numbers of blades as shown in Table 


(I). Therefore, 
jones, 5. 4, Seeeoewor. | (3522) 


However, the two bladed propeller, i.e., 2 = 2, is not very 
common in conventional merchant and naval ship designs and, 
therefore, is not included in this study. 

2) Equivalent Reynolds Number--The Series data, as pub- 
lished, is valid only in the range of equivalent Reynolds 


numbers given by: 


2x10 eae! Zeit Coa) 





If the equivalent t/c is varied from the standard equivalent 


value (t/c) 75R! then the new equivalent Reynolds number 


i ee) which results from this variation, must lie within 


the same limits. That is, 
2 <0 


2x10 (see) 


5 


aOR 


These limits are appropriate for full-size propellers. For 
example, given the following: 

a) wt = .22 

bye = 20 (knots) 

C) Np = 104 (rpms) 

eo = 25. (ft) 


e) c arr 0 (Ube) 


.75R 
£) v = 1.2285x10 (£t*/sec) 


a 


the value for Rn fsseguaieiae) 2.0.9 x 10”. 


Be go) 


3) Pitch-diameter Ratio (P/D)--The series data considers 


only pitch diameter ratios in the range given by: 
Oo, P/D =< 1.4 (3425) 
4) Advance Ratio (J)--An inspection of the Series results 
in graphical format shows that J varies over a range given 


ey < 


OR ecmegar< ING (3.26) 
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5) Expanded Area Ratio (A,/A,) -~Using Table (I), AL/Ag 
varies over certain ranges depending on Z. This is stated 


aS: 


0.35 < A,/A, < 0.8 Z = 3 (S27) 
0.40 < A,/A, < 1.0 Z= 4 ie 28) 
0.45 < Ap/AQ < 1-05 Z= 5 onc) 
0.50 < A,/A, < 0.8 Z = 6 (3330)) 
oe a =) 08S ig Gael) 


6) Hub diameter-to-Propeller Diameter Ratio (d/D)--From 
Table 37, Section 17 of Reference [16], the Series data 


requires that: 


apy = O16 Z 


iI 
WW 
~ 
~~] 


(34-32) 


O26 7 Z 


il 
i 
~ 
oR 
~ 
>) 


(es) 


i, 


a/D 


foe 2 UMMARY 
From the preceding discussions, the following observa- 
tions can be made: 
1) the "Design Cases", defined in Chapter I, are examples 


of the powering equation solution approaches. That is, Design 
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Case No. 1 constitutes a "Thrust Approach" problem; Design 
Case No. 2, a "Power Approach" one; Design Case No. 3, a 
"Matching" problem. 

2) equations (3.17) and (3.8) imply that an optimization 
solution to the "Design Cases" will involve Poe ews eel 
P/D, Ap/Ag: Cele n: Qq and Np as possible design variables. 

3) when viewed from the concepts on optimization presented 
Piechapter II, equations (3.25) through (3.31) constitute 
Side constraints to an optimized solution of a propeller 
selection problem which uses the Wageningen B-Screw Series. 
Having noted these points, the propeller selection problem 
can now be formulated as a general, non-linear, constrained 


optimization problem. 
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ie enor nsth ee onotBeL LON=-AN OPTIMIZATION PROBLEM 


A. INTRODUCTION 

The purpose of this chapter is to present the formula- 
tion of the propeller selection problem as an optimization 
problem that can be solved using COPES/CONMIN. Three usual 
restrictions considered by the designer in any propeller 
selection analysis are stated as constraints. Then, the 
components of D and X are assembled based on requirements 
from previously cited relationships. The restrictions 
considered by the designer and the limitations imposed by 
the use of the Wageningen B-Screw Series are presented in 
inequality constraint format. A formal statement of the 
propeller selection problem as an optimization problem is 
followed by a review of the GLOBCM common block format and 
the basic subprograms used in all three versions of SUBROUTINE 
ANALIZ that pertain to each Design Case. 

The FORTRAN subprogram listings are found in Appendix 
B. Comment cards have been used extensively in the coding 


development to assist the reader. 


B. DESIGNER'S CONSIDERATIONS 


ieenopel ler 512¢ 
The first restriction on the selection of any 


propeller is size. That is, the propeller race in the stern 


of the hull under consideration will only accommodate a 


a2 





propeller of some given maximum diameter (Diam): As a con- 
straint on a selected propeller of diameter D, this may be 


written as: 


un 
or, alternatively, as: 


oe = 0 (4.2) 


Sa 





2. Cavitation 

Another item of importance in propeller selection is 
the cavitation phenomenon. When a propeller of given diameter 
D and expanded area ratio Ap/Ap 1S operating to produce a 
thrust T, the formation and subsequent collapse of water vapor 
bubbles on the blade surface, i.e., cavitation, is likely to 
occur if the localized surface pressures, usually on the 
"back" side of the blade, drop below the pressure at which 


the fluid would boil (p ) in the surrounding environment. 


watvap 
Avoidance of cavitation can be reasonably assured by selecting 
a propeller having certain geometric characteristics. A good 
empirical relationship that establishes these characteristics 
for propellers typified by the Wageningen B-Screw Series is 
the Keller Cavitation criterion [Ref. 2: p. 259]. It speci- 


fies the minimum required expanded area ratio (AL/A co 


Oe tion 


avoid cavitation and is given by: 


Bis 





Choa ae, 37) - 
(AJA)... = ee — tb (4.3) 
eget (Patm tT Pacd hay Pwatvap’ D* 


where: 
Z = number of blades; 
aD = developed thrust (l1bf); 
D = propeller diameter (ft); 
Pe em = atmospheric pressure (psia); 
Pwatvap Twin vApOrizazLomepressure (psia) 
2 4 
QO = fluid density (lbf sec /ft ) 
2 
acg — Ze) 4 {ftysec ) 
Dag acre Mmro shart centerline (£t) 
b —SeCOustant-s Ua! tor 2 =.2, 0.2 for 2 = 1. 


As a constraint on the propeller selection, this requirement 


1S writen as: 


(AE/AG) min Ap/Ag oo 
Or 
CR PRAY oh 
= _ Oech oe 


3. Strength 
The final designer's consideration (for this study), 


maeluded in the selection of a propeller, is that of strength. 
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Given the propeller's material (promat), selected from Table 
(II), and the loadings (T and Q.) imposed, it is important to 
ensure that the blade's cross sections have proper dimensions 
(in an ideal sense, maximum blade section thickness (t*) and 
chord length (c)) to ensure adequate strength. Since the 

use of the B-Screw Series requires that the chord length 

(c) vary aS a prescribed function of D, Z and Ap/A,: as given 
in Table 1 of Reference [2], the adequacy for strength can 

be determined by an appropriately selected value for blade 
section maximum thickness-to-chord ratio (t*/c) alone. So, 
slag Ee is the established minimum blade section maximum 
thickness, then the strength requirement follows from the 


constraint given by: 


eas 
min 


= * * 
Ve Ge [a ese cc (4.6) 
The fact that blade section maximum thickness for the B-Screw 
Series varies linearly with the propeller radius (R) allows 
the strength constraint (equation (4.6)) to be evaluated at 
one section along the radius. This point is chosen to be at 


the 3/4 radius (.75R). Therefore, equation (4.6) becomes: 


ve*/°) 75R min < ‘t°/°) 5p Serot) 


Or 


(t*/¢) 75R min 
aa CD 


aD 





Reference [2] suggests the following empirical relation for , 


the minimum required equivalent blade section maximum 


thickness~-to-chord ratio (t*/c) 75R ea 





(Ze) Seo, Dy) 


D 


ZecOmOOes+O.2 0 


2 ee 
4.123n,D°(s, +2 Np 
.) a eZ oS 
aI, a POW mene A 
: oe) 
(4.9) 
where: 
D = propeller diameter (ft); 
Py = delivered power (hp); 
N, = propeller revolution rate (rpm); 
- = propeller material allowable stress (psi); 
P/D = pitch-diameter ratio. 


However, in Chapter VI of this thesis, an algorithm which 
employs the Schoenherr formulation [Ref. 8] with some modi- 


fications, is presented as an alternative to equation (4.9). 


fee ink DESIGN VECTOR 

In view of the preceding presentations on optimization 
and powering, the design vector D can be assembled for the 
general propeller selection problem utilizing the B-Screw 
Series. This vector is composed of preassigned parameters 


relating to environmental conditions, hull characteristics 
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and the propeller which are required for various equations 
and the design variables. 
1. Parameters 
a. Environmental 
These parameters pertain primarily to the fluid 
conditions in which the propeller operates and to the atmos- 
phere. Required for various calculations, they are: 
1) fluid temperature (°F) --Temp 
2) fluid density (lbf sec*/ft")--p 
3) fluid viscosity (ft*/sec)--v 
4) fluid vaporization pressure (psia)--p 


watvap 
5) atmospheric pressure (psia)--p_., 
b. Hull Characteristics 
These parameters pertain to certain details pre- 
scribed for the hull under study in the powering analysis. 
They are: 
1) wake fraction--wt 
2) thrust deduction--td 
3) relative rotative efficiency tl 
4) number of propellers--noscrw 
5) shaft centerline depth (ft)--h_, 
6) propeller diameter limit (ft)--D, 
ere rope. ler 
These parameters are specified in view of their 
discrete-valued nature. They are: 


Ll) number of blades--Z 


2) material--promat 


=) 4) 


Table (II) lists materials and properties considered in this 
study. These values are taken from Table (35), Section 15 


of Reference [16]. 


TABLE II 


Material Identifier Reference 


promat Material Allowable Stress--Sc Density--wd 
ee (1lb£/in?) 
i Cash iron SoU 0—— ae 0 . 260 
2 Cast Steel 5915-~-6265 ig 
3 Type 2 Bronze 7200-7585 Po 
4 Type 4 Ni-Al 8910--9430 Pea 
Bronze 
5 Stainless Steel 5400--5500 w233 


eeebesign Variables 
ineeyvrew Of Equations (3.13) through (3.17), the 


design variables common to all selection approaches are: 


i) P, 

may 6OUV 
3a) 

4) P/D 
5) An/AG 


6) (t*/c) sep 
7) N 
8) Q 


The vectors D and X are shown schematically in figure (4.1). 
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Figure 4.1 


Design Vectors D and 


5g 
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D. CONSTRAINTS 

Besides the constraints imposed by equations (4.2), 
(4.5) and (4.8), equations (3.25) through (3.31) are rearranged 
to the format of constraints in equation (2.9). They are 
listed as follows: 


1) Equivalent Reynolds Number--Equation (3.25) becomes: 


* 


Rn 
1 =—=—— = 1000 (4.10) 
2x AO 
Two constraints are derived: 
* 
Rn 
G,(X) = 1- —23 0 WL) 
2 GeO 
ee 
6,(X%) = —+22- i000 < 0 (4.12) 
2x LO 


2) Expanded Area Ratio--Equations (3.27) through (3.31) 


become: 


(AB/AQ) lower '7) = Ap/Ap = (AE/Ad) upper (2) 
(A513) 
Two constraints are derived: 
Ge (X) = (AE/A,) lower (2) - Ap/Ag <0 (4514) 
Ge (X) ~ AL/AG = (AE/A,) upper ‘2? < 0 (4.15) 
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3) Advance Ratio--Equation (3.22) becomes: 


0 < =—- < l cao) 


Ree eee 2 
Cue ee = eC (Ane) 
2 Ta x 


4) Equivalent Blade Section Maximum Thickness-to-Chord 


Ratio--Using equation (3.19), boundaries on the range of (t*/c) I5R 


are defined by: 


B(t/o) asp < (t*/e) g5q < 41t/c) ap dae! 


Two constraints are derived: 


LA 
© 


ze l 
a 7S) Gor” «('t 7S) a5R eo 


1A 
>: 


Gp (X) = (eee (an 1) 


.75R ~ *4t/C) a5R 
ObJECTIVE FUNCTIONS 

Upon consideration of equation (3.13), Design Case No. l 
and Design Case No. 2 require that the open water efficiency 
ae) given by equation (3.11), be maximized. In terminology 


related to optimization, this is stated as: 
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OBJ, 2 rea) (4.22) 
Design Case No. 3, the "matching" problem, requires that the 


blade weight (bldwt) be minimized. This is stated as: 


OBJ = bldwt (4.23) 


F. PROPELLER SELECTION OPTIMIZATION PROBLEM STATEMENT 
As a general, non-linear constrained optimization problem 
to be solved by COPES/CONMIN, the propeller selection problem 


for all Design Cases may be stated as one equation given by: 


Minimize: F(X) = OBJ, , or OBJ, (4.24) 
Subject to: G(X) 0 i= lap atet lle 
xe =) ee 1 6 
aL — i- 1 
mae constraint Gy 5 (X) and the values for Keener and ea 


will be specified according to each Design Case. 


G. CODING FUNDAMENTALS 
1. GLOBCM Common Block 
The GLOBCM common block, required by COPES/CONMIN, 
is now assembled. Table (III) specifies the assignment loca- 
tions for the FORTRAN variables which define objective 
functions, design variables and constraints. 
2. SUBROUTINE ANALIZ 
While each Design Case uses a different approach, 
all analyses are very similar. Therefore, each SUBROUTINE 
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TABLE IIT 


Global Common (GLOBCM) Catalog 


Global FORTRAN DEFINITION 
Location Name 
1 ETAO is 
2 WEIGHT bldwt 
3 AEDVAO AD/Ay 
4 DIA D 
5 N Np = 60 ny 
6 PE Py 
7 PDIVD P/D 
8 QS on 
* 

9 TC75R (t*/C) sep 
10 V Wet 7 Sec) 
ll RJCHL G, (X) --eqn Aa) 
2 RICNU G, (X) --eqn (4.18) 
13 R75RCL G(X) --eqn TA) 
14 R75RCU Cede ecdn G42) 
15 AEAOCL G, (X) --egn eae) 
16 AEAOCU G, (X)=-eqn (4215) 
17 Meno CL G. (X) --eqn (4.20) 
18 TC75CU Gp (X) --eqn 2 )) 
19 POWBAL G,.(X)--eqn (7.10) or (8.11) 

eZ 

or (9.4) 

20 DIACNU Gg (X) --egn (4D) or (9). 6) 
21 AEAOCV G19 (X) --eqn 2S) 
22 TCSTRS G, 1 (X) --eqn (4.8) 
23 Rel J 





ANALIZ shares a common structure and other common subroutines 
which perform calculations required in all cases. Appendices 
C, F and I contain, respectively, the source ieeeas or 
SUBROUTINE ANALIZ for Design Case No. 1, DeSign Case No. 2 
and Design Case No. 3. 
oe oO eLEuCcCEure 
The structure common to all cases follows 

pi@cordingly: 

ieeeall initialization of environmental, hull and propeller 
parameters is accomplished in the input section (ICALC = l). 

2) evaluation of Kan and Ko: 


objective functions (—n 0 or bldwt) are accomplished in the 


all constraints and appropriate 


execution section (ICALC = 2). 
3) output of results for each optimization problem is 
accomplished in the output section (ICALC = 3). 
b. Basic Subprograms 
The following FORTRAN subprograms are used in 
all three SUBROUTINE ANALIZ codes: 

1) SUBROUTINE CH75RA--calculates the equivalent blade 
section chord length (C sep) for the propeller using Table 1 
mor. 2, p. 252). 

2) SUBROUTINE REY75R--calculates the equivalent Reynolds 
number (Rn*~.5) using equations (3.19) and (3.20). 

3) SUBROUTINE COEFSA--calculates the thrust and torque 
coefficients (Ky, and Ko) through sequential calls to SUBROU- 


TINE CALCKT and SUBROUTINE CALCKQ. The polynomial expressions 
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(Tables (5) and (6), [Ref. 2]) for these coefficients are 
contained in SUBROUTINE CALCKT and SUBROUTINE CALCKQ respec- 
tively. . 

4) SUBROUTINE OPWEFF--calculates the open water efficiency 
(no) Mising equation (3.11). 

5) SUBROUTINE JCNA--calculates the constraints on the 
advance ratio (J) given by equations (4.17) and (4.18). 

6) SUBROUTINE REYCNA--calculates the equivalent Reynolds 
number constraints given by equations (4.11) and (4.12). 

7) SUBROUTINE EXTCCN--calculates the constraints on ex- 
panded area ratio (AE/A)) and equivalent blade section maxi- 
mum thickness-to-chord ratio (t*/c s.p) given by equations 
Meera), (4.15), (4.20) and (4.21). 

8) SUBROUTINE DICNUA--calculates the constraint on the 
propeller diameter (D) given by equation (4.2) using the 
hull parameter on maximum diameter (Dy. ) > 

9) SUBROUTINE CAVCNA--calculates the constraint for 
cavitation given by equation (4.5) using equation (4.3). 

10) SUBROUTINE STRCNA--calculates the constraint for 


strength given by equation (4.8) using equation (4.9). 


H. SUMMARY 

The propeller selection problem has now been formulated 
aS a constrained optimization problem which can be solved 
by COPES/CONMIN. Two items remain for discussion before 
proceeding to specify the final details pertaining to each 
SUBROUTINE ANALIZ code and to present numerical examples. 


These items are: 
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1) the theory and coding relating to the computation of 
the propeller's blade weight for the evaluation of the objec- 
tive function in Design Case No. 3 (OBJ 5) . 

2) the theory and coding relating to the computation of 
the minimum required equivalent blade section maximum thickness- 


to-chord ratio (t*/c 75 a) for use in the alternative 


evaluation of the strength constraint given by equation 


(4.8). 
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V. PROPELLER BLADE WEIGHT--AN OBJECTIVE FUNCTION 


A. INTRODUCTION 

In this chapter, the method for the computation of the 
propeller's blade weight (bldwt), the objective function 
OBJ, is examined. First, a brief overview on the steps in 
the computational procedure is presented. The FORTRAN 
subprogram SUBROUTINE WGTCAL developed from the algorithm is 
then described. Again, Appendix B contains all subprogram 


Pstings. 


B. THEORY AND PROCEDURE 

Given the material of a propeller, the calculation of 
the weight of one blade involves nothing more than a volume 
calculation, a relatively routine task performed by most naval 
architects/marine engineers. Analagous to the determination 
SE, che underwater volume of a ship's hull, the calculation 
1s an integration of blade section profiles' cross-sectional 
areas over the propeller radius (R). 

Perum es Of Integration 

Figure (5.1) depicts a side elevation view of a 

blade and hub, parallel to the propeller shaft axis. The 
cross-hatched area indicates the trace of the volume to be 
Calculated. In view of equations (3.32) and (3.33), limits 
@aeincegqration are from r = .167R to r = R for Z = 4,5,6 and 
r= .18R to r = R for Z = 3,7. For convenience, a non- 


dimensional variable "x" will be defined as: 
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Figure 5.1 Propeller Blade & Hub--Side View 





(5) 


AIK 


MiG TOrewtoste x = 1.0. 


Limits are now expressed as x 

Quite obviously, R= (D/2.0). 
2. Blade Section Profile 

Figure (5.2) depicts an expanded cylindrical blade 

section in profile view at a given r or x. For the Wageningen 

B-Screw Series, the profile is defined, geometrically, by 

a succession of vertical ordinates which specify points along 

the blade section's profile on the "face" (Y¢) and on the 

back (y,) with respect to the pitch reference line. At any 

r = xR, vertical ordinates for "aft" (P < 0) and "fwd" 


= 


(P > 0) portions of the blade section are determined by: 


eed = Nabe ™ > te.) Bane (Se2) 

TO Cime(tta.) +t, P< 0 (E23) 
and 

ee Vy (t* = fe Eee 0 (5.4) 

ee Gy v5) (E* — t*) ae Pen (5.5) 
where: 

Vi +V>5 = tabulated values depending on x and P 


(see Tables (2) and (3), Ref. [2]); 
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* 


vo - blade section leading edge thickness (ft); 
ae ode msoCemOlmenailang cdge thickness (ft) . 


Marts for Yea Yoa’ Yee and Yp¢e are feet (ft). For this 


Study, a reasonable assumption is made in that: 
_ E IL * 
1S = t = (+a) © (5:76) 
3. Blade Section Cross-Sectional Area 


The cross-sectional area at each x = r/R (A(x)) is 


determined by: 


9 10 
BASS) ae) evan (Se?) 
i=l i=l 
where 
ioe +h, 
_ ai+l al 
Dog ES 3 (5.8) 
a ~ bai ~ Yai eee 
et ee me eee page eee 
Expressions for Aey Nes and Desay follow in similar fashion. 


Values for Yeni and Yoai ore determined at 9 points 
along the "aft" portion of a given blade section's chord (c): 


values of Yeri and Ypeq are determined at 10 points along the 


al 





"fwd" portion. The values for AP; and AP es are fractional 
values of the blade section's chord length (c) at radius 

r = xR as determined from Tables (2) and (3) in Rerorence 
[2]. The units for A(x) are square feet (ft*). Units for 
h h are feet (ft). 


mh c, AP_. and AP 
al 


ai ot” «fa — £itl’ ie 


4. Volume Integration 
The blade volume (bldvol) is finally determined by 
using Simpson's Rule for integration of A(x) along the non- 
dimensional radius x uSing appropriate limits. 
5. Blade Weight 
Once the blade volume (bldvol) is calculated, the 


weight (bldwt) is determined by: - 
bldwt = bldvol - wd . 1728 Coe etal.) 


where: 


bldvol = volume of one blade fo 


wd material weight density (Ib£/in”). 


Weight Density (wd) depends on blade material (promat). 


Table (II) lists appropriate values. 


ee CODING 
SUBROUTINE WGTCAL is the main subprogram for the blade 
weight calculation. It, in turn, calls the following FORTRAN 


subprogram for various calculations: 
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1) SUBROUTINE TDIST--generates, at specified radius 
values, a distribution of blade section maximum thicknesses 
mes). . 

2) SUBROUTINE BLDPRP--generates, at specified radius values 

Suet — .lo7R or .18R, .2R, .3R, .4R, ..., .9R, 1.0R), 
various "blade section properties", one of which is a blade 
section's cross-sectional area given by equation (5.7). 
Other properties which are determined (for later use in 
direct stress computations) include blade section chord 
lengths and centroids and "critical point" locations as 
defined in Chapter VI. 

3) SUBROUTINE BLDVOL --performs a Simpson's Rule integra- 
tion for the propeller blade volume (bldvol) using blade 
section cross-sectional areas generated in SUBROUTINE BLDPRP. 
The blade weight (bldwt) is computed as a final step in the 
main subprogram SUBROUTINE WGTCAL. 

Examination of the codes in Appendix B reveals extensive 
use of common blocks for passing data from one subprogram to 
another. Comment cards provide a full definition of all 
common blocks as well as a description of the task being 


performed at various points in a given subprogram. 


D. SUMMARY 
The coding developed for this study is, admittedly, not 
very compact and efficient. However, the intention has been 
to write all codes with sufficient documentation in order to 
facilitate the reader's understanding of the algorithms employed 


as well as to make the author's debugging work easier. 
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Viera ene oo-GO—-CHORD RATEO--A DESIGN CONSTRAINT 


ee) 6«LNTRODUCTION 

The purpose of this chapter is to examine the development 
of an algorithm that will be used to determine the minimum 
required equivalent blade section maximum thickness-to-chord 
ratio used in equation (4.8). The formulation is based upon 
the method developed by Dr. Karl E. Schoenherr [Ref. 8] in 
1963. After a review of the past and present methods employed 
in propeller strength analysis is conducted, a description of 
the Schoenherr model and a list of the assumptions used with 
that model is presented. Then, a brief restatement of his 
model's equations which are used in the algorithm is followed 
byea derivation of the author's modifications to the Schoenherr 
method. The chapter is completed by conducting a review of 
the theory and coding employed by the algorithm. 

The principal reference which is cited throughout this 
chapter is, again, Reference [8]. The reader is encouraged 


to review this reference for further details. 


B. PROPELLER STRENGTH ANALYSIS--A HISTORICAL REVIEW 

Marine propeller blades present a special class of struc- 
tural problem. That problem lies in the difficulty of des- 
Cribing a blade design in simple mathematical terms for 
subsequent analysis through various means. Until the "finite 


element method era", analytical methods, including the one 
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[Ref. 8] adapted for this study, relied heavily on practical 
experience of the propeller designer and semi-theoretical 
considerations. Analysis by these methods provided a cri- 
terion of stress rather than actual computation of stresses. 
These methods for predicting blade stresses were developed by 
using "beam" theory or "shell" theory. 

The use of elementary beam theory in propeller strength 
analysis was first adopted by Taylor [Ref. 22]. He treated 
a blade as a cantilever beam attached to the propeller hub 
and loaded by thrust and torque forces distributed linearly 
over the propeller radius. His approach is often deemed a 
"modified beam theory" because he chose to calculate the 
direct stresses using the moment of inertia properties of 
expanded cylindrical blade sections with neutral axis parallel 
to the nose-tail (pitch-reference) line or chord of that 
expanded section. Reasonable estimates of stresses along the 
blade surface were achieved for the unraked, unskewed and 
narrow-bladed propellers of his time. 

As propellers "modernized" and became skewed and wider 
with increasing rake (mostly aft), modifications, improvements 
and alternatives to Taylor's theory were developed. Prin- 
Cipally, modifications by Rosingh [Ref. 23] and Hancock 
[Ref. 24] proposed using moment of inertia properties of a 
blade section that was normal to the generating line of the 
axially projected blade outline. Romson [Ref. 25] later 


improved Taylor's theory for application to wide-bladed 
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propellers. Morgan [Ref. 26] provided an improved method 

for calculating the geometric properties of “modern" airfoil- 
shaped blade sections. Aernoldus and Keyser's [Ref. 27] 
"“guasi-static" modeling of the propeller blade allowed 

for additional consideration to stresses induced by centri- 
fugal loading of the raked and skewed blade. The beam 

theory approaches to propeller blade stress analysis culminated, 
for all practical purposes, with Schoenherr's work [Ref. 8] 

ol 963. 

Alternatives to Taylor's beam theory approach, prior to 
1963, consisted of the application of "shell" theory to the 
propeller blade strength problem. This approach was first 
proposed by Conn [Ref. 28] and subsequently formulated by 
Cohen [Ref. 29] who modeled the blade as a helicoidal shell 
with variable thickness and infinite width. "Shell Theory" 
was utilized again in experimental studies by Connolly 
[Ref. 30] who, like his predecessors, was also forced into 
making an assumption about the behavior of the displacements 
of the blade sections (i.e., constant displacements normal 
to the constant pitch blade at each fractional radius dis- 
tance from the hub) beyond usual assumptions of shell theory. 
Essentially, his experimental results on one specific pro- 
peller contradicted the computational values. Attempts at 
a generalized numerical solution to Connolly's equations 
appeared in 1963 [Ref. 31] and 1964 [Ref. 32]. In 1968, 


Atkinson [Ref. 33], compared Connolly's results with currently 
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adopted cantilever beam methods and found, based on the 
inconsistency of results, that it was not possible to. 
recommend one method over the other in the blade strength 
design procedure; another approach was needed. 

Commenting on Atkinson's paper at that time, Sontvedt 
[Ref. 34] pointed out that, in view of these inconsistencies, 
the only approach to blade strength analysis which did not 
require very broad assumptions was the finite element tech- 
nique. Developments in the method at that time were provid- 
ing a new and powerful tool for structural analysis. The 
propeller blade was just another application. Genalis [Ref. 
35] developed codes for the determination of displacements 
and stresses in a blade under hydrodynamic loads using the 
FEM technique and modeling the blade as a shell, a 3-D ele- 
ment mesh of tetrahedrons and rectangular prisms and, finally, 
a composite of shell and 3-D elements. As an aside, a 
finite "difference" solution to Connolly's analytical equa- 
tions was proposed in 1972 [Ref. 36]. In 1973, Atkinson 
[Ref. 37] reported the application of both hydrodynamic and 
centrifugal loads to a blade modeled by a thin-shell triangu- 
lar mesh and a thick-shell parabolic and cubic curved element 
mesh. The results of the triangular element were considered 
unsatisfactory. Another use of the thin-shell triangular 
element was reported by Sontvedt [Ref. 34] in 1974 using the 
SESAM-69 code [Ref. 38]. 

The need to model the blade correctly near the hub, 


where root stresses are usually critical, necessitated the 
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consideration of 3-D elements in lieu of thick/thin-~-shell 
elements [Refs. 39,40]. The use of the 4- and 10~noded 
tetrahedral element (i.e., TET~4, TET-10 respectively) to 
construct a blade mesh was conducted by Beek [Refs. 41,42]. 
He observed that improved accuracy of stress values, achieved 
by the use of these meshes, were overshadowed somewhat by 

the extensive storage capacity required for each analysis. 
Another “natural" improvement, from a geometrical standpoint, 
appeared in 1978. A general 3~D curved isoparametric element 
was incorporated in a computer code [Ref. 43] developed by 

Ma based on his previous formulation work [Ref. 44]. 

-. The finite element approach will continue to grow in use 
in propeller blade strength analysis with each successive 
improvement made to the basic elements which are used in the 
mesh generation of the blade. But, until the computer storage 
problem is resolved to the point where one mesh generation 
and subsequent stress analysis of one particular blade be-~ 
comes a minor processing task, the basic analytical techniques 
will continue to be a meritable "check" [Ref. 45] in the 
preliminary (or conceptual) phase of propeller selection/ 
design. In this context, the method formulated by Schoenherr 
and his colleagues twenty years ago is considered for adoption 
in determining a minimum required blade section maximum 


thickness~-to-chord ratio. 
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C. SCHOENHERR'S METHOD 
ie background 

Schoenherr's method is applicable to meee 
(or, conceptual) propeller selection problems because it 
employs an assumed thrust and torque force loading distribu- 
tion for the propeller blade. This assumption is made at 
this stage of the ship's design because the exact wake 
velocity distribution at the propeller race is generally 
not known. 

Also, his method is applicable to propellers repre- 
sented by the Wageningen B-Screw Series because the blades of 
Emese propellers meet Schoenherr's criteria for the blade 
types covered by his formulation. Specifically, B-Screw 
Series blades have: 

1) a small constant rake angle of 15° over the entire 
blade radius; 

2) a constant pitch distribution over the propeller 
radius with the exception of the Z = 4 propeller whose pitch 
is slightly reduced near the hub; 

3) mild skew; 

4) linear distribution of blade section maximum thickness 
over the radius of the blade; 

5) aerfoil profile qualities where the nose-tail line 
or the chord of the blade section is approximately parallel 


to the pitch reference line. 


Ve, 





2. The Blade Model 
Selecthere models) the propeller blade as a cantilever 
beam with unsymmetrical and variable area cross sections 
subjected to loading distributions of hydrodynamic and cen- 
trifugal forces. The following additional assumptions apply: 

1) Flexure theory applies. This subsequently implies the 
following: a) plane cross sections remain plane under load, 
b) Hooke's Law is valid, c) the blade material is homogeneous 
and isotropic, d) fibers are free to extend and contract 
independently of adjacent fibers, and e) stresses at a point 
arising from various forces superimpose. 

2) Shearing stresses and their effects are neglected. 

Only the direct stresses on a strength section are taken 
mace account. 

3) The strength sections are taken to be the expanded 
cylindrical blade sections at various radial locations. 

4) The neutral axes of a strength section are straight 
lines passing through the centroid of the expanded cylindrical 
blade section and are parallel and normal, respectively, to 
the pitch reference line, and therefore, the chord, at each 
blade section. 

5) Bending Moments are applied in two planes which are 
mutually perpendicular to each other. One plane is normal 
to the pitch-reference line (and chord) of the strength 
section; the other is parallel. 

6) The angle between the principal axes of inertia and 


the neutral axes is zero. 
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Using this model and assumptions, Schoenherr applies 


the following formula for the evaluation of the direct fiber 


stress eo) in a blade section at a radius r = ro: 
[M] u [M] W ee) 
a ne © von Eo 
lol, = i + r + Gan (6.1) 
no Q 
where: 
eo! IM} og = resultant bending moments in planes 
normal ([M] a and parallel CIM]. EO 
the strength section's chord at r= ry 
(Get ahon a; 
a! 4 = centrifugal force acting normal to the 


plane of the strength section at r=r 
and resulting from the centrifugal 
acceleration of the remaining blade 
element mass above that strength section 
Cig); 


a, WwW = coordinates of a point on the strength 
section's periphery with respect to 
that section's neutral axes system (2-n 
system) (ft); 


A(x_) = strength section's cross-sectional area 
. (EES) 


I = moment of inertia of the strength 
section with respect to the "2" axis 
ete, 


a = moment of inertia of the strength section 
WiheMegespect co the “n™ axis (f£t*) >; 


xe = non-dimensional radius given by 


as r/R: 


Since equation (6.1) indicates that the direct fiber stress 
1S greatest at points on the periphery of the strength sec- 


tion, Schoenherr selects to examine four "critical points" 
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on the periphery where the fiber stress is likely to be a 


maximum. These points are designated as (2), (2), (3) and 


(4) on Figure (6.1) and are specified by coordinates (ul,wl), 
(u2,w2), (u3,w3) and (u4,w4) respectively in the "=n" 
reference system. 


The values for ae and [M]} are determined by 


LO 


the following relation: 


[M] = [M_] + [M 


no P=ne BG (6.2) 


[M] = [(M_] + [M 


ce P° Lo oes! 


EyeiiS 


where: 


= total bending moment due to hydrodynamic 
loading acting in a plane normal to a 
strength section's chord at r=r 
(ftlbf) ; ° 


Dene 


[Mp1], = total bending moment due to hydrodynamic 
O ' : 
loading acting in a plane parallel to 
a strength section's chord at r = ry 
(ie lo) =. 


= total bending moment due to centrifugal 
loading acting ina plane normal to a 
strength section's chord at r=r 
(Coe Neva) ae 


BAG 


O 


BS 06 = total bending moment due to centrifugal 
loading acting in a plane parallel to 
a strength section's chord at r=r 
(£tlb£) - : 


3. Bending Moments Due to Hydrodynamic Loading 


The derivations of ene and [MoJo, EOulow GdLlrect ly 


from Part I of Schoenherr's paper [Ref. 8: p. 83-89] and, 
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therefore, only key equations will be restated. References 
to equations which contain no decimal point apply to equations 
as numbered in his paper. 

Thrust and torque force are components of the hydro- 
dynamic "lifting" force acting on a blade. Using an assumed 
non-linear distribution of thrust along the blade radius 


given by equation (2), Schoenherr derives the following ex- 


pression for the bending moment due to thrust ((M, J.) which 
acts at a blade section located at ratius r = ro: 
6,(x_) 
7 2 
a. eee en 
Sah 
where: 
7 = propeller thrust (lbf); 
R = propeller radius (ft); 
Z = no. of blades; 
d5 (x ) >, (%)) = functions of non-dimensional 
: radius x evaluated at x_ = r/R 
and x, = 0.2 and given By equa- 
tilows 4) ana. (9). 
For the bending moment due to torque (IM 1) which 
acts at a blade section located at radius r = los Schoenherr 


derives the following: 


Q Ue) 
P 2 iO 
[M_] = —_—_— +s (6. 5} 


where: 
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Q Seon pelten EQrque (ft—-1bf) ; 
z, = MQ, 2O ir blades; 


Jets), o.(s ) = functions of non-dimensional radius x 
2 ae. Pn 
evaluated at x, =r /R and x, = 0.2 
and given by equations (es?) Qnd Cay). 


Figure (6.2) depicts the component resolution for 
[Mp1] 50 and IModog which results when equations (6.4) and 
(6.5) are imposed at a strength section at r = ry which has 
pepitch angle Bo: The following relations are derived as 


equations (42) and (43) in Schoenherr's paper: 


(M,] = (M 


De eels cos a + IMJ, sin EP (6.6) 


IMalog = : cos B (6.7) 


[M ea Sain B - Mody 


4. Force and Bending Moments Due to Centrifugal Loading 


The derivation and expressions contained in this 
section constitute the author's modifications to the formu- 
tation in Part II of Schoenherr's paper. In Part II, Schoen- 


herr's derivations for [M and [M are formulated for 


eee 


computation uSing a propeller drawing. This follows from 


Bone 


the fact that his method, which was funded by the American 
Bureau of Shipping, was intended to be used as that classi- 
Fication society's "designer's check" on adherence to the 
Bureau's strength criteria from a propeller blueprint. To 
evaluate the direct fiber stresses from equation (6.1) for 


the Wageningen B-Screw Series in accordance with Schoenherr's 
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] 


method, [M and [M must be evaluated from expres- 


Gb: 20 


mens derived from the available information on blade sec- 


Een 


tion profiles and other geometric characteristics which 


are contained in Reference [2] and previously used in Chapter 


V. 

Consider Figure (6.3) where the centrifugal force 
[Ch], of a blade element above a blade section at x = x, = ey 
acts in a radial direction from the shaft centerline. Its 


line of action passes through point "N", which is on the same 
cylindrical surface as the blade section at x = aS r/R, 
and through point "G", which is the center of gravity of the 
blade element above the blade section at x = eae r/R. 


From the figure, the follcwing expression is derived: 


Flo = [C1], cos es) tc [Cp], sin ce) ere) 


]. 1s shifted to point "N" and is decomposed into 


[C,, z 


components [C.] cos ae) and [Ch] sin C(x). The entire 


O 
cylindrical surface in which the blade section at x = —. r/R 
and point "N" lie is now expanded into a flat plane for further 
consideration (see Figure (6.4)). In this configuration, 
(Cal, cos ¢(x,) is normal to this flat plane while [C,], sin ¢(x,) 
lies in this plane. 

Let point "O" be the location of the blade section's 
neutral axes system (i.e., the &£-n system). Then, the forces 


and moments due to the centrifugal reaction of the blade ele- 


ment above this section act at point "O" and are given by: 
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eos = [C1 cos g(x) (6.9) 
PDedG = [C,], sin é (x) (65 10) 
eee s Lae - NO (Geet) 
[Mowe - Oe i yO (6212) 
where: 
ree — sGitece, terce, Gue Lo Centrifugal action, 
S acting on the blade section located at 
X = x = x _/R; 
O O 
De! 5 = shear force, due to centrifugal action, 
acting on the blade section located 
at x = x_ = xr _/R; 
O O 
Ee 6 = bending moment at point "O" imposed by 
Piel Ace rMiGacIrOoulghn point 'N"; 
re = torsional moment at point "O" imposed by 


Lele Ageing cIrougiepOint 'N'. 


Since Schoenherr's method does not consider shear forces and 


their effects, [D,] and [M_.] will not be considered in 
O cw’o 


this modification. However, eye. and eee must now be 
computed for each blade section along the propeller's radius 
in order to account for their contributions to equations 


eee, (6.3) and, finally, in equation (6.1). 


The computation is derived as follows. Consider 
Figure (6.4). Again, Eel acts through point "N" and is 
normal (outward) to the plane of the figure. Ua lis is now 


resolved into components of the 2-n axes system as follows: 
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[Mplno = [Cpl cos 5(x,) tp, sinBo +q,cosB +¥ Qo} (6.13) 
[Mop] oo [cn], cos o(x,) {q. sin BOW Po cos BG + mes coer 4) 
where: 
ae = pitch angle of the blade section at 
Koso = © o/ Re 
O O 
Bee = distance of the blade sections centroid from 
the generator line (ft); 
a = distance of the blade section's centroid from 
the pitch reference line (ft); 
= dastamec c@ cOornt "N” from point "Pp" 
parallel to the shaft axis at x = xo = r/R 
(ES ; 
De Gustance to point "N” from point "PP" 


perpendicular to the shaft axis at 
a eS r/R (Git). 


The quantity Bg TSmeouncd oy the relation: 


She 


tan B = (P/D) |, (6.4.15) 
For the Wageningen B-Screw Series, 72) is a constant along R 
except for propellers with Z = 4. 


The quantity [C IP 1s computed from the relation: 


EF 
wd V D _ 
me) = 1728 ° ee Sn) ae (XC -) (6.16) 


where: 


oul 





V = volume of the blade element above the blade 
section at x = x = r/R (ft3); 


acg = 32.174 ft/sec’; 
wd = material weight density (lbf/in?) ; 
n = propeller revolution rate (rps); 
(X,) = non-dimensional radial position of "G" 
: for the blade element above the shaft 
axis; 
R = propeller radius (ft). 


At this point, only five quantities remain to be determined 
for the evaluation of the expressions of equations (6.10), 
(6.13) and (6.14). They are: 


1) cos ¢(x,) 


2) By 
3) qo 
4) (%) 
5) vs 


These quantities are determined by integration over the 
blade element above the blade section, located at x = ae r/R, 
from x = x Ome = 2,0. 

The values for Py and qo will vary with the location 
of "G" (from Figure (6.3)) which depends on x. Consider a 
radially thin slice of the blade element above the blade 
section, located at x = x, = r/R (see Figure (6.5)). This 
thin "slice" is located at a non-dimensional distance x from 
the shaft centerline where se 1.0. Figure (6.6) depicts 


this section expanded onto a plane. Let "g" be the centroid 


Seeeaat “slice”. I£ a 1s the distance from the generator 


oz 
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imme FO "g" in feet and ee is the distance from the pitch 
reference line in feet, then, from Figures (6.5) and (6.6), 


the following relations apply: 


5 = x R tan n Corel?) 
1 -a = 1 -x_ sin - cos 6.18 
g g g g Pg 7 Pg \ 
ie = xX_ cos = Sin Gspacikies 
g g Pg  ¥q Pg 
where: 
nN = rake angle at x; 
oe = distance of "g" from point "P" parallel to 
J the shaft axis (ft); 
y = distance of "g" from point "P" perpendicular 
J to the shaft axis (ft). 


For the Wageningen B-Screw Series, rake angle n is a constant 
15° everywhere along the radius R. 

Now, to compute the volume of the blade element above 
the blade section located at x = Sa r/R, the following 


expression is used: 


V = J RA(x) dx (6.20) 


To compute the non-dimensional radial position of 
meertor the blade element above a blade section located at 


a. = r/R, the following expression is used: 
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J 


f Rec ieo dx 
X=X 2 
ec) = ——____ (ce2) 
Go 1 ; 
f A(x) dx 
x=X 
O 
To compute the tangential position of "G" for the 
blade element above a blade section located at x = x_ = r/R, 
the following expression is used: 
1 
i A(x) t_ dx 
XSX_ 
TS = T G22) 
ii A(x) dx 
X=X 


And, finally, to compute the axial position of "G" 





for the blade element above a blade section located at x = Xo 
- r/R, the following expression is used: 
il 
A(x) (i_-a_) dx 
J (x) (ig-a, 
eae 
A, T KO. 3) 
f A(x) dx 
x 
Oo 
Using the values just determined for (Xa) o T, and 
A, the following expressions are used to evaluate Py and 
do ° 
*o 
De = == at (6.24) 
ees) : 
Go 
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aan A - XS tan n Gi) 
The expression for cos G(X,) follows: 
cos G(X) = —— (6226) 


The formulation is now complete. Equations (6.10), 
(6.13) and (6.14) can now be evaluated for any blade section 
located at x = a r/R. From here, equations (6.2) and 
(6.3) are evaluated. Finally, using the results from these 
equations and equation (6.9), equation (6.1) can be evaluated 
for the four points, specified by coordinates (ul,wl), 
(u2,w2), (u3,w3) and (u4,w4), at any location x = B= r/R. 


From the development discussed in Chapter V, the 


Sees for A(x), x_, y., I and I are readily available. 
g g no LO 


D. ALGORITHM FOR THE CONSTRAINT 
1. Theory 
Schoenherr's formulation with the modifications just 

derived can be used in determining the minimum required 
equivalent blade section maximum thickness-to-chord ratio 

# e . Ts * 
ect /C) s6p min) for use in the constraint Gj, (%) < 0 given 
by equation (4.8). The procedure employed is as follows: 

Soe ie ; 

1) assume an initial value for (t /C) sep min Using 

equation (3.21); 


2) increase this value by a small amount; 


oF 


Be using (t*/c) obtained from step (2), generate 


7 ee 
a distribution of minimum required blade section thicknesses 
(tran? for blade sections at specified points along the 
Smepeller radius, say at r = .2R, .3R, .4R, .5R, .6R, .7R, 
.8R and .9R; 

4) determine all blade section properties to include: 
a) cross-sectional area, b) chord length, c) centroid location, 
d) moments of inertia with respect to the principal axes 
system (1.e., £-n system), e) coordinate values for the four 
critical points defined in the previous section; 

5) compute the hydrodynamic bending moment components 
and [M 


[M at radius locations just specified; 


Eno p! go 
6) compute the values of the centrifugal force epee and 


= 


Eales acting on 


the bending moment components [M and nie 


oe 
blade sections at radius locations just specified; 
7) calculate the direct fiber stresses at all four 
critical points for all radius locations specified in step 
(3); 
6) check the following condition on the calculated fiber 


@eeess at all four critical points at all specified radius 


Mecations using: 


io < 144 + § fo 27) 


9) if the maximum allowable stress (S.) for the material 
is exceeded, then return to step (2) and repeat. Otherwise, 


proceed to next step. 
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10) since the minimum required equivalent blade section 
maximum thickness-to-chord ratio assumed in step (2) has 
produced blade sections of adequate strength, evaluate the 
constraint given by equation (4.8). 

2. Coding Details 
The algorithm just outlined is incorporated into the 
main FORTRAN subprogram SUBROUTINE STRCNK. This subprogran, 
in turn, executes the algorithm through sequential calls to 
other key FORTRAN subprograms. These subprograms are listed 
as follows: 

1) SUBROUTINE TDIST--accomplishes step (3); generates, 
at specified radius values, a distribution of minimum required 
blade section maximum thicknesses (tr an) for the assumed 
value of (t*/c) yep win) 

2) SUBROUTINE BLDPRP--accomplishes step (4); described 
previously in Chapter V; 

3) SUBROUTINE HYDLD--accomplishes step (5); computes the 
hydrodynamic bending moment components, given by equations 
(6.6) and (6.7), at specified radius locations; 

4) SUBROUTINE CNFGLD--accomplishes step (6); computes the 
centrifugal force and bending moments, given by equations 
(6.9), (6.13) and (6.14) respectively, at specified radius 
Hecations; 

5) SUBROUTINE SIGNDS--accomplishes step (7); computes 


direct fiber stresses, given by equation (6.1), for all four 


Critical points at every specified radius location. 


hs, 





During the remaining steps of SUBROUTINE STRCNK, the 
condition on allowable stress, given by equation (6.27), 
is checked at all critical points of blade sections located 
at specified radius locations (again, r = .2R, .3R, .4R, .5R, 
.6R, .7R, .8R and .9R). The final calculation made is 
that for the constraint given by equation (4.8). 

Again, extensive use of common blocks, for passing 
data from one subprogram to another, is apparent upon examina- 
tion of the codes just cited. Comment cards are used 


Peroughout. 


E. SUMMARY 

The end of this chapter marks the completion of all 
prerequisite background and formulation discussions on the 
application of COPES/CONMIN to propeller selection problems 
involving the Wageningen B-Screw Series. From this point, 
each specific Design Case can now be solved as an optimiza- 


tion problem. 
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VII. DESIGN CASE NO. 1--PROGRAMMING AND COMPARISONS 


foe LNTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the “thrust” approach. 
First, the thrust approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solution 
to this problem is presented. Four variations to this pro- 
peller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 


moe results from the four variations. 


B. THRUST APPROACH FORMULATION 
1. Design Vector Xl 
AS previously pointed out at the conclusion of Chap- 
mere.lil, Design Case No. 1 constitutes a propeller selection 
problem which is solved by the thrust approach. In this ap- 
proach, the effective horsepower (Pp) and the ship's speed 
(V) are specified by the designer. From the viewpoint of 


optimization, the quantities P., and V become preassigned 


2 

parameters. This reduces the design vector X (see Figure 4.1) 
Lo : 

D 

P/D 

xX = AD/AG (Zak) 
* 
ese 
Np 


OEE 





Having specified P._, and V, all of the design variables, 


EB 


as listed in equation (7.1), are not independent. Recalling 


Semeacions (3.3) and (3.13), the following relationship 





results: 
Pp = {itstd) = ° N o anes o a) (7 2) 
E (l-wt) OR O 5510 60 ' 
Rearranging terms, this equation becomes: 
feo 550160 
Reo 2 N See OE Cis) 


o “S 'P (1-td) 27 


Ue: 
Considering that the open water efficiency one) is evaluated 
prior to the computation of aD; Ors OBJ; 9: then both N5 
and Q,. are not independent design variables. One must 

be selected as the independent design variable. Then, 

the other variable becomes dependent on the one just 
selected. 


For this study, N, is selected as the independent 


Ie 
design variable. This choice will reduce the design vector 


X for propeller selection problems using the thrust approach 


to the following: 


P/D 
X = AL/A ied) 


NaC SR 
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Finally, equation (3.8) implies an alternative defini- 
tion of X as given in equation (7.4). The design vector for 
Design Case No. 1 propeller selection problems is, therefore, 


defined as: 


P/D 
XI = ( A,/A, (7.5) 
SSeS 


2. Powering Constraint 


Having determined the design vector Xl, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 
restriction constitutes the remaining constraint Gy 5 (X) men- 
tioned in Chapter IV. 

Simply stated, the selected propeller, as defined by 
X1, must develop enough thrust (T) so that the powering require- 


ment, specified by P, and V, is met. Using equation (3.2), 


E 
the thrust developed by the propeller can be specified in terms 


of thrust horsepower (Pin) as: 


T V(l-wt) 


(Pm) Goy = 55 (iG) 
fmeemeequation (3.9), it follows that: 
0 ne D4 Kn 
ydesis 550), © “°) oe 
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Moeng equation (3.12), the developed thrust horsepower can 


Bemderined in terms of developed effective horsepower given 


Dy : 
Pec). 
Pa? sen ~  (l-wt) (Pn) dev ao) 
The restriction imposed by the thrust approach method, 
where P.,, and V are specified, can now be stated as: 


E 


: = ( 


ace (PS) (7.9) 


dev 


Rearranging equation (7.9), the constraint Gj 5 (XI) follows: 


Cap) a ee C76) 


With the design vector X1 and G, (XI) defined, the 
propeller selection problem represented by Design Case No. l 


can be stated under one equation as: 


Maximize: F(X1l) = OBJ 


Sua ject to: Grea =) 70 Ae ei ies yt 2 (oy aed) 


Pee REVLOUS SOLUTIONS 
Triantafyllou [Refs. 3,21] considered a propeller selection 
problem represented by Design Case No. 1. In his example 


problem, the following parameters were specified: 
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2 (eee cee 


Meese < 108° (m-/sec) =0%.22613 x 107 
2) Wie = 2 2 
Mied = .19 


4) O22 


NR 


5) noscrw = lL 


6) 22 5 

i) Pe = 18153 (hp) 

8) V= 24 (knots) 

fee D = 22.0 (£t) 
10) Ap/Ay aS 


The hull under study in his example had the following dimen- 
sions: 
ime Length = 710 (ft) 
meme Draft = 30 (ft) 
3) Beam = 100 (ft) 
For his analysis, the design vector contained two varia- 


bles and was specified as: 


P/D 


Using an iterative scheme [Ref. 21: p. 79] to solve two equa- 
tions in two unknowns, he maximized the open water efficiency 


=5! to obtain the following results: 


P/D ae 65-1 


N 


p 104 (rpm) 


EOS 





aw) 
il 


25544 (hp) 


.6676 


3 
ll 


For future comparisons, equations (3.8) and (3.3) give: 


.8286 


J 


Qe eee Cr Ome feast ) 


Triantafyllou's results are summarized in Table (IV). 


D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
meet), 1S now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 

The first and second variations attempt to reproduce the 
solution given by Triantafyllou. The design vector XT 
(NDV = 2) 1s used in both cases. One variation uses SUBROUTINE 
STRCNA to evaluate the constraint G5 (XT) given by equation 


(4.8). The other uses SUBROUTINE STRCNK to determine G,.(XT). 


12 
The remaining two variations will solve the propeller 
selection problem using the design vector Xl (NDV = 5) defined 
in equation (7.5). Again, one variation uses SUBROUTINE 
Bene NA; the other, SUBROUTINE STRCNK. 
In all variations, the following parameters are used: 
meetemp = 59 (°F) 
2) 9 = 1.9384 (lb£-sec*/£t*) 


Mme 1.2285 10 (£t/sec) 
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4) Pwatvap = .247 (psia) 
5) oa 14.7 (psia) 
oS wt = .22 


om td = .19 


O23 


8 ) Np 


9) noscrw = l 


nO) h =e OF EE) 
CY 

Pl) Po am eno) (it) 

me) 4= 5 


13) promat = 5 (stainless steel; see Table (II)) 

14) Pa = ess (ap) 

moe V = 24.0 (knots) 
All of the above are initialized in the input phase (ICALC = 1) 
of each SUBROUTINE ANALIZ pertaining to each variation. 

The constraint Gy 5 (X1) or Gy 9 (XT) is evaluated by SUBROU- 
TINE BLPOW1L which appears in the execution section of each 
SUBROUTINE ANALIZ. 

ieee Variation l 

a. Programming Details 
SimGemtibism variation uses the design vector XT, 

- the following design variables of Xl become parameters and 
are specified in the input section of SUBROUTINE ANALIZ (ICALC = 1) 
as: 

ieeeD = 22.0 (ft) 

Zz) A,/AG = .85 

Bb) (t*/c) se, = .0348 (from equation 3.21). 


Oe 





For constraints, the following are used: 
G, (XT) < 0 i= lee 8, 12 


Only nine of twelve constraints are evaluated (NCON = 9). 
Obviously, some of the twelve constraints are redundant since 


b, A,/A poand (t*/c} have been specified. 


POR 
upper 
a 


lower 


) and lower (XT ances On 


The upper (XT 


the design variables J and P/D are set to be: 


ao - + .o 
wes | by < 1.4 


These upper and lower limits are specified in the COPES 
control card deck on card image F under respective fields 
VYUB and VLB. The initial value for each design variable 
(XT, ) is also assigned on card image F under the field labeled 
feeeine first list of card images in Appendix D lists all of 
the COPES control cards used for this variation and variation 
2. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this varia- 
tion, found in Appendix C, shows the calling statement made 


EeeoUSROUTINE STRCNA. 
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b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed first in Appendix BE. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 
ye Variation 2 
a. Programming Details 
Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix C, shows 
the calling statement made to SUBROUTINE STRCNK. 
b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed second in Appendix E. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (IV). 
fee Variation 3 
a. Programming Details 
This variation uses the design vector Xl. For 


constraints, the following are used: 
G. (x1) Os) Gab 35 5 2 


All twelve constraints are evaluated (NCON = 12). 


The upper Oa ld cies and lower (xj tower 
ve cs 


) ima. Es 
. ; 2 
on the deSign variables D, P/D, Ap/Ao: ie /C) sep and J 


are set as: 


O'S 





ieee P< 50.0 (ft) 
eee 7 Oe 1 4 : 
-2 < APA, < tel 
* 
Ce /Cle mes 50 


ee <0 


These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (Xl, ) 

is also assigned on card image F under the field labeled xX. 
The second list of card images in Appendix D lists all of 

the COPES control cards used for this variation and variation 
4. These cards also specify the locations of the design 
variables in the common block GLOBCM (see Table (III)) as 
well as the locations of the constraints and their boundaries. 
Further details on the COPES control card requirements and 
the format of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix C, shows the calling statement 
Made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed third in Appendix E. Re- 
sults for this variation of thepropeller selection problem 


are tabulated in Table (IV). 
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4. Variation 4 
a. Programming Details 
Everything discussed above for the third varia-~ 
tion applies here with one exception. An examination of 
SUBROUTINE ANALIZ for the fourth variation, found in Appendix 
C, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA. 
b. Results 
The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed last in Appendix E. Re- 
sults for this variation of the propeller selection problem 


are tabulated in Table (IV). 


iO LSCUSSION 

Overall, the results achieved in all variations compare 
reasonably well to the solution obtained by Triantafyllou. 
However, the following points can be made. 

Variations 1 and 2 give the same results. This was ex- 
pected in view of the fact that, even though constraints 
Gy (XT) throuch G, 4 (XT) were evaluated, these constraints were 
not considered in the optimization search conducted by CONMIN. 

In variations 3 and 4, the diameter (D) was driven to 
the limit (Dian): This bears out a fundamental rule in pro- 
peller design, i.e., the larger the propeller diameter (D), 
the greater the open water efficiency (No) - 

The minimum required equivalent blade section maximum 


thickness-to-chord ratio ((t*/c) jeeecompukea in 


.7/5Rmin 


eae 


variation 3, is substantially smaller than the one computed 
for variation 4. As pointed out in Reference [2], the empiri- 
cal relation, expressed by equation (4.9) and derived from 
equation (70.x) [Ref. 46: p. 620], does not take into account 
the effects of centrifugal loading. These effects include, 
specifically, the direct stresses imposed by the inertia load 
of the blade and the bending moments which result from rake 
and skew of the blade. Therefore, the algorithm developed in 
Chapter VI should, and does, produce a larger value for 
7") 7S5R° 

A final observation on the results concerns the values 
of the open water efficiency. The "optimum" open water effi- 
clency (no) achieved by Triantafyllou is lower than those 
achieved in variations 1 and 2. A possible reason for this 
might be the neglection of the term "dRe/dJ" in Triantafyllou's 
Peemulation of the analytical expressions [Ref. 21: p. 71) 
that he used in his analysis. The difference in the open 
water efficiencies subsequently accounts for the differences 
in the propeller revolution rate (N.,) and the delivered 


torque (Q.) when the relation in equation (3.3) is considered. 


gee 





Design Case No. 


GROUP ITEM TRIANTA- 
Pyne) 
Po 1352.0 
Given V Zao 
Design D Za.0 
Variable 
Bpeci— AR Ao a 
; * 
fied (Ge /C) 357 
1D) 
Design P/D iGo: 
a A/a, 
* 
ne7<) oR 
J . 8286 
Maximize ae .6676 
Da im 
Restric- 
tions comin 
. 75Rmin 
N, 104 
Other Q. 1290000 
Py 755440 


TABLE IV 


1l--Results 
1 2 
£353 .0 6153.0 
24.0 24.0 
22.0 22 0) 
85 aBo 
.0348 .0348 
1.0036 00 36 
Bs A. esl 
— OSL ~/091 
(2250) (2250) 
(0.5256) (0.52538) 
(0221266) 0 250701) 
116.9 iG. 9 
1080451. 1080451 
Z405i.3 24051 .3 


Hales 


5 


Neila -0 
24.0 


ZA 9991. 
Beheiok 
2205 
70.3350 
. 7343 


a LO2 


Zon 
p00 
sOZ923 


117.4 
1064574 
240 30 .3 


4 


oes 0 
24.0 


PEN ENS yo)s, 
iO 7 
8149 
0642 
. 1394 


- /066 


2280 
poco 
nUSozoU 


ioe / 
1084003 
24094 .8 





VIII. DESIGN CASE NO. 2--PROGRAMMING AND COMPARISONS 


mee) «= LNTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
propeller selection problems which use the "power" approach. 
First, the power approach to the propeller selection problem 
is formulated. Then, a review of a previous author's solu- 
tion to this problem is presented. Four variations to this 
propeller selection problem are solved by COPES/CONMIN. The 
chapter is completed with a presentation and discussion of 


the results from the four variations. 


B. POWER APPROACH FORMULATION 


1. Design Vector oO) 


AS previously pointed out at the conclusion of 
Chapter III, Design Case No. 2 constitutes a propeller selec- 
tion problem which is solved by the power approach. In this 
approach, the delivered torque (Q.) and the propeller revolu- 


tion rate (N.) are specified by the designer. From the view- 


ie 


point of optimization, the quantities Q. and N. become pre- 


E 


asSigned parameters. This reduces the design vector X (see 


meoure (4.1)) to: 


< 


X = D (8.1) 
P/D 
An/Ag 
oR 
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Having specified Q, and N all of the design varia- 


Ep 
bles, as listed in equation (8.1), are not independent. Re- 
Belling equations (3.3), (3.11) and (3.13), the following 


relationship results: 





cee eS) Poe: (8.2) 

E (I-wt) ° "R 27K, © 550 60 
Rearranging terms, this equation becomes: 

i - ii-td) , , fi-wt) . au esi Np Gene) 

V (I-wt)  "R “nl D K 550 0 

1 Q 

Considering the relations for Kon and Ko imewecuatwton (3217) ; 
then both P_ and V are not independent design variables. 


E 


One must be selected as independent, while the other becomes 
dependent on the one selected. 

For this study, V is selected as the independent de- 
Sign variable. This choice will reduce the design vector X 
for propeller selection problems using the power approach to 


the following: 


V 

D 

P/D (8.4) 
Ap/AG 

(t*/C) a5p' 


e 
—— 


<I 
I 


Finally, equation (3.8) implies an alternative defini- 


£ion of X as given in equation (8.4). The design vector for 


LS 





Design Case No. 2 propeller selection problems is, therefore, 


defined as: : 


V 
J 

XD = (| P/D ioe 
Ap/Ao 
(E*/C) o5R 


2. Powering Constraint 


Having determined the design vector X2, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 


restriction constitutes the remaining constraint G,.(X) men- 


iy 
tioned in Chapter IV. 

Simply stated, the selected propeller, as defined by 
X2, must absorb at least all of the power delivered to it 
(Pp) which is specified in terms of Qq and N,- Using equation 
(3.3), the power absorbed by the propeller can be specified 


in terms of delivered horsepower (P)) as: 





Zee N 
(P_) ~ ee (8.6) 
D‘ absorb 550 6 : 
Beomeequation (3.10), it follows that: 
K. 0 me D> 27 N 
(p_) oe eee (8.7) 
D‘ absorb 550 60 ° 


But, equation (3.3) also defines the power delivered to the 


propeller as: 


ln 





auc 


D 550 °&#&«&©60 





8re8) 


The restriction imposed by the power approach method, 


where Qq and Np are specified, can now be stated as: 


P < 


D — Sy Aer (8.9) 


Rearranging equation (8.9), the constraint Gy 5 (X2) follows: 


(P_) 
pee —Psesoxe Zea (8.10) 
D 


meeener Simplification of equation (8.10) gives: 
Gj 5 (X2) = b= — <- 0 Corde) 
With the design vector X2 and G15 (X2) defined, the 


propeller selection problem represented by Design Case No. 2 


can be stated under one equation as: 


Minimize: F(X2) = OBJ, , 
Subject to: G5 (X2) SO Ga elgg) 2 (cee) 
ee eee ee t= 1,...,5 


fee REVLIOUS SOLUTIONS 
Markussen [Ref. 4] considered a propeller selection prob- 


lem represented by Design Case No. 2. In his example problem, 


x] 





the following parameters were specified: 


1) Temp = 18 (°C) = 64.4 (°F) 


Z) Pwatvap = 0.0206411 (bars) = .29943921 (psia) 
3) a = 1.01312856 (bars) = 14.6974 (psia) 

4) noscrw = l 

>) Ao = 6.7 (meters) = 21.9827 (ft) 

6) Z2= 6 

7) Pp = 18.9 (MegaWatts) = 25344.9 (hp) 

8) Np = 110 (rpm) 

feey. = 15.65 (knots) 


For his analysis, the design vector contained three 


variables and was specified as: 


a | 
P/D 


Ap/Ag \ 


2 


A restriction for the minimum required expanded area ratio 


((A,/A ), given by equation (4.3), was also considered. 


o’ min 
This imposed a constraint given by equation (4.4). 
Using an iterative scheme [Ref. 4: p. 110] to solve three 


equations in three unknowns, Markussen maximized the open 


water efficiency (n) to obtain the following results: 


J = .61095 
P/D = .864380 

ei 2 
Ap/Ag =oGeelool/40.6123 (m /m ) 


io 





Sree bz 


~o04 39:1 


No 


meretuture comparisons, equations (3.8) and (3.3) give: 


O 
I 


feeeuotetmeters), = 23.593375 (ft) 


IL JL HOI) (nga olan) 


tO 
I 


Markussen's results are summarized in Table (V). 


D. SOLUTIONS BY COPES/CONMIN 

The propeller selection problem, as stated by equation 
(8.12), is now solved by COPES/CONMIN. Four solution varia- 
tions are considered. 

The first and second variations attempt to reproduce the 
Beegcion given by Markussen. The design vector XM (NDV = 3) 
is used in both cases. One variation uses SUBROUTINE STRCNA 
to evaluate the constraint Gy 5 (XM) given by equation (4.8). 


The other uses SUBROUTINE STRCNK to determine G (XM) . 


12 

The remaining two variations will solve the propeller 
selection problem using the design vector X2 (NDV = 5) de- 
fined in equation (8.5). Again, one variation uses SUBROUTINE 
STRCNA; the other, SUBROUTINE STRCNK. 


In all variations, the following parameters are used: 


1) Temp = 64.4 (°F) 


mmo = 1.9892 (1bf-sec*/ft') 
fe = 1.1900 x10°° eee eee) 
4) Pwatvap = .2994 (psia) 


iD 





5) = 14.697 (psia) ! 


Patm 
6) wt = .22 


me tad = .19 


a ON £2025 


R 
9) noscrw = lL 


0) Aa —e Zi oo2 / (ft) 
il ) em = SHU SC ag) 
ma 24 = 6 


13) promat = 5 (stainless steel; see Table (II)) 


Poor sOnfe— 1p t ) 


14) Qo 


i>) N JEJE) aaron, 


i 
All of the above are initialized in the input phase 

(ICALC = 1) of each SUBROUTINE ANALIZ pertaining to each 
variation. 

The constraint Gy 5 (X2) or G, 5 (XM) is evaluated by 
SUBROUTINE BLPOW2 which appears in the execution section of 
each SUBROUTINE ANALIZ. 

1. Variation l 
a. Programming Details 

Since this variation uses the design vector XM, 


the following design variables of X2 become parameters and 


are specified in the input section of SUBROUTINE ANALIZ 


(ICALC = 1). The ship's speed (V) is specified as: 
V = V,/ (1-wt) 
—Meeeomooy (lL = ,22} 


20.0641 (knots) 


EZ 0 





Markussen elected to use the standard Wageningen blade sec- 
tion maximum thickness-to-chord ratios. Since the equivalent 
t/c is given as a function of Z and AB/AG (see equation (3.21)), 


Saenm (t*/c) 55R 1s calculated during each analysis (ICALC = 2) 


by the following relation: 


CES/C) a5p = (t/C) sep 


For constraints, the following are used: 


G, (XT) = i eles, oy. Ores 12 


Only ten of twelve constraints are considered (NCON = 10Q). 


Constraints Gy (XM) engl (@ 1 (XM) are redundant since no limit 


i: 


Pameene propeller diameter (D ) appears aS a parameter in 


lim 
Markussen's formulation and (t*/c) 75p was taken to be the 
Wageningen standard. 


lower 


oe.) and lower (XM" 


The upper (XM 5 \_ ELM ES 


on the design variables J, P/D and Ap/A5 are set to be: 


rol =S J eee 


4 S P/D 


[A 
rae 
na 


A S Ap/AG <P bat 


These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 


and VIB. The initial value for each design variable (XM. ) is 


eek 





also assigned on card image F under the field labeled xX. 

Gage first list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 2. 
These cards also specify the locations of the design varia- 
bles in the common block GLOBCM (see Table (III)) as well as 
the locations of the constraints and their boundaries. Further 
details on the COPES control card requirements and the format 
of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, 
performed by COPES/CONMIN, is listed first in Appendix H. 
Results for this variation of the propeller selection problem 
are tabulated in Table (V). 

fee Variation 2 
a. Programming Details 

Everything discussed above for the first variation 
applies here with one exception. An examination of SUBROUTINE 
ANALIZ for the second variation, found in Appendix F, shows 
the calling statement made to SUBROUTINE STRCNK instead of 
SUBROUTINE STRCNA. 

b. Results 
The output from the optimization/analysis, per- 


formed by COPES/CONMIN, is listed second in Appendix H. 


ee 





Results for this variation of the propeller selection problem 
are tabulated in Table (V). 
3. Variation 3 
a. Programming Details 
Piemvanmiatienmeises the desian vector X2. For 
constraints, the following are used: 


G; (X2) < 0 fea) ee 2 


All twelve constraints are evaluated (NCON = 12). 
The upper (Age) and lower one.) limits 
on the design variables V, P/D, AD/AG: tee) I5R! and J are 


set as: 


Om Os V < 100.0 (ft/sec) 
Ss P/D i 1.4 
4 < AL/A5 Ss ila! 

MOU ee ce) (E“7C) oop < ou 
5 Wl J < eee 


These upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X2.) 

is also assigned on card image F under the field labeled xX. 
The second list of card images in Appendix G lists all of the 
COPES control cards used for this variation and variation 4. 


These cards also specify the locations of the design variables 


2S 





in the common block GLOBCM (see Table (III)) as well as the 
locations of the constraints and their boundaries. Further 
details on the COPES control card requirements and the format 
of each card are contained in Reference [7]. 

An examination of SUBROUTINE ANALIZ for this 
variation, found in Appendix F, shows the calling statement 
made to SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed third in Appendix H. Re- 
sults for this variation of the propeller selection problem 
are tabulated in Table (V). 

4. Variation 4 
a. Programming Details 

Everything discussed above for the third varia- 
tion applies here with one exception. An examination of SUB- 
ROUTINE ANALIZ for the fourth variation, found in Appendix 
F, shows the calling statement made to SUBROUTINE STRCNK 
instead of SUBROUTINE STRCNA. 

b. Results 

The output from the optimization/analysis, per- 
formed by COPES/CONMIN, is listed last in Appendix H. Results 
for this variation of the propeller selection problem are 


tabulated in Table (V). 


fee OLoCUSSION 
The results achieved in variations 1 and 2 compare ex- 


tremely well to the solution obtained by Markussen. As 
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Pemiteca GUL in the discussion in Chapter VII, variations l 
and 2 are expected to give the same results for the vector 
XM. Obviously, the values obtained for J and P/D, as well 

* 
menthnose for D, and Rn sspe 


generated in Markussen's example. However, the values for 


are very close to the values 


AL/Ag are somewhat different. It is interesting to note that 
the value obtained in variations l and 2 (and, for that matter, 


variations 3 and 4) is, essentially, the limiting value for 


AB/AG: as given in Table (I), for Z = 6. Markussen's value 
for AD/AG (i.e., .891012) exceeds the limit (i.e., .80) in 
this table. 


As pointed out at the end of Chapter VII, the minimum 
required blade section maximum thickness-to-chord ratio 


)) , computed in variations 1 and 3, is sub- 


me 7°) 75R min 
Stantially smaller than the one computed for variations 2 
and 4. Again, the same explanation applies here as well. 
The results of variations 3 and 4 differ somewhat from 
Markussen's results. The reason for this is simply that Vis 
(or V) has not been specified as a parameter. Consequently, 


a higher value for the advance ratio (J), which corresponds 


to a higher open water efficiency (No) s has been found in 


the optimization search. This result can be interpreted 
in the following way. Given: 

1) a six-bladed Wageningen propeller (Z = 6) which is made 
Out of stainless steel (promat = 5); 


2) a power train delivering 25344.9 (hp) at a rate of 


EO (rpm) ; 


ie S 





3) a hull with a wake fraction (wt) equal to .22, a 
thrust deduction (td) equal to .19 and a shaft centerline 
depth (h_,) of eo ee), . 
then, the selected propeller, as defined by X2, can drive 


this hull at a maximum speed of V when the hull has a maximum 


resistance given by Poe 


6 





TABLE V 





Design Case No. 2--Results 
GROUP ITEM MARKUSSEN VARIATIONS 
1 2 2 4 
P, 25344.9 25344.9 25344.9 25344.9 25344.9 
Given oF ieee oNOnsO 1210030 1210130 1210130 
N,, ielG 110 110 110 110 
Design 
Variable Va 15.65 15.65 15.65 
Scecified 
a V DOmiGdl  S20m0G4 1 
a 61095 6475 6475 9927 8753 
V 30.9138 29.5355 
Design V. 2a 27. Z3anOore. 
Variables 
P/D 864 380 9036 9036 1.1986 1.0308 
AA, 891012 8018 8018 7946 7986 
* 
(ene) aan 0397 0397 0499 0638 
Maximize Ny 654391 6660 Facce 7643 7330 
oe (30.0) (30.0) 30.0 30.0 
Restrictions A,/A,_., 574729 5070 5070 4622 vas 
* 
(e-7/e) oe (.02729)  (.0647) 02706 0638 
D Da 52 D225 ee Sm 2286 24.23 
Other Pe (14057.3) (14057.3) 19945.5 20653.7 
Rn* 3 6.478x10 5x10/ es 10) 8 es calo 
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IX. DESIGN CASE NO. 3--PROGRAMMING AND COMPARISONS 


eae LNTRODUCTION 

In this chapter, COPES/CONMIN is used in the solution of 
a propeller selection problem where "matching" is desired. 
First, the "matching" approach to the propeller selection 
problem is formulated. Then, a review of a previous author's 
solution is presented. One variation to this propeller 
selection problem is solved by COPES/CONMIN. The chapter is 


completed with a presentation and discussion of the results. 


B. "MATCHING" FORMULATION 
1. Design Vector x3 
Design Case No. 3, the final powering problem con- 
Sidered in this study, constitutes a propeller selection prob- 
lem solved by the "matching" approach. In this approach, 
the hull's effective horsepower (P.) and speed (V), the 


delivered torque (Q.) and the propeller revolution rate 


S 


(N,) are specified by the designer. This reduces the design 


vector X (see Figure (4.1)) to: 


D 
P/D 

Apso 
eer 


| 
i 


(a) 


For this study, the design vector X is reduced further 


by eliminating the propeller diameter (D) as a design variable. 


Ws 





That is, D will also be specified by the designer so that 


fe design vector for Design Case No. 3 is defined as: 


P/D 
oa AQ/AG (9.2) 
NES) aia) 


Z. Powering Conseraint (s) 


Having determined the design vector X3, a final 
restriction to the general propeller selection problem, as 
stated by equation (4.24), remains for consideration. This 


restriction constitutes the remaining constraint G,.(X) men- 


~ 


le 


tioned in Chapter IV as well as an additional constraint. 

In the "matching" problem, the selected propeller, 
as defined by X3, must satisfy two conditions. First, it 
iaee develop, aS a minimum, the effective horsepower (P.) 
as imposed by the design specification. Citing the formula- 
tion previously derived in Chapter VII, this condition can 


be stated as: 


: PR = (PE) dev (9.3) 
The constraint G15 (3) follows accordingly as: 
(PS) 
Cee 7 oe! eG (9.4) 
2 PS = 


For the second condition, the selected propeller can 


only absorb, as a maximum, the delivered power (Pp) as 


ano 





specified by the designer. The formulation is the same as 
that in Chapter VIII except that the inequality signs are 


reversed. The condition is stated as: 


E (97 5)) 


Ee Abeer = D 


moederining a constraint G13(X3), qneeie pelocat ton, say 

location 24, in the GLOBCM block (see Table (III)) would be 
assigned. But, considering the fact that constraint Gy (X) 
will not be used because the propeller diameter (D) is speci- 


fied, there is no reason why Gy (X3) Cannot be redefined, for 


this Design Case only, as: 


(P..) 
(X3) = __D absorb _ 1 an} (9.6) 
Py _ 


Gg 


Further simplification of equation (9.6) gives: 


(X3) = ~-1 < 0 Cee) 


In reality, the constraints just defined should be 
equality constraints. The word "match" does infer equality 
in some sense. However, as previously stated in Chapter II, 
the version of COPES/CONMIN used in this study does not 
directly handle equality constraints. But, since CONMIN 
attempts to minimize constraints in the optimization search, 


meme! be assumed that a "match" can be achieved. 
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With the design vector X3 and the constraints 
G15 (X3) and Gy (X3) defined, the propeller selection problem 
represented by Design Case No. 3 can now be stated under 


one equation as: 


Minimize: F(X3) = OBJ, 
Subject to: G, (X3) < 0 ape pss aise 2 (9.8) 
x3lower <a < ahae a a= a lie no) 
aa = = ae 


eRe VIOUS SOLUTIONS 

The propeller selection problem considered by Vassilopoulos 
[Ref. 18} actually represents a propeller "design" problem 
using the "power" approach. The example problem which he 
elected to solve is taken from that posed by the International 
Towing Tank Conference (ITTC) Propeller Committee. This 
problem is concerned with the determination of propeller 
thrust (T), diameter (D) and speed of advance (V.) (or ship 
speed (V)) for a single-screw cargo ship where: 


1) power available to the propeller (i.e., Py) is 


e0,000 (hp) 
2) Z2= 6 
EyeeN = 105--110 (rpm) 
5) Nog = 19 (ft) 


The variation of ship speed (V) and of hull effective power 
(Po), thrust deduction factor (td) and the wake fraction 


(wt) is also given [Ref. 18: p. 20}. 
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The results from Vassilopoulos'’ propeller design exercise 
Beoauced a propeller that is "matched" at the following values: 
i) Pa = aie 2.0 (hp) 
mee Vv = 24.24 (knots) 


3) Qo = S00 606.70 (LE—l Dr) 
4) Np = 105 (rpm) 
5) P) = 30000.0 (hp) 


His propeller "design" was based upon the following specified 
Parameters: 
im Temp = 59 (°F) 


Mm = 1.9905 (1bf-sec/£t*) 


3) Pyatvap = .247 (psia) 
meewe = .22 

mee ca = .1725 

6) noscrw = l 

7) Ao = 19;0 (£t) 

8) Z= 6 


9) promat = 5 (stainless teel, see Table (II)) 
foe = 22.0 (ft) 


el) Np 


105 (rpm) 


i) Pp SUCOOs 05 (hp) 


D 
Using an optimization scheme incorporated in his MVAPDP 
computer program, Vassilopoulos maximized the open water 


efficiency ee and designed a propeller with the following 


Characteristics: 


SZ 





3) K, = .0478 


Q 
4) No = .,691 
5) Ap/Ap = .767 
mm Dladwt = 7617.2 (1bf) 


By utilizing both the lifting line and lifting surface 
methods in his design procedure, Vassilopoulos' MVAPDP program 
Seyorvyed a "Constant stress” propeller blade. Consequently, 
the values for (t*/c) and P/D varied non-linearly along the 
propeller radius (R). According to Vassilopoulos, this 
resulted in a minimum weight propeller. The values for P/D 
and (t*/c) are listed in Tables (8) and (10) of his paper. 
From these values, (t*/C) s5p 1S approximately .040. 

While the propeller represented by Vassilopoulos' desigr 
is different, in many aspects (rake, skew, blade section 
aerfoil shape, etc.), from the Wageningen B-Screw Series 
propeller, it does represent a minimum weight propeller that 
has been "matched" to specific design values. Appropriate 


results are summarized in Table (VI). 


Pee SOLUTIONS BY COPES/CONMIN 
The propeller selection problem, as stated by equation 
(9.8), is now solved by COPES/CONMIN. One solution variation 
is considered. The following parameters are used: 
meee cemp = 59 (°F) 


Ics Sm ibesssc- 7c) 


2) p= 

=> 2 
Sey = 1.2817 x 10 Clee 77sec} 
4) Pwatvap ~ -247 (psia) 


esi 





5) Eee, 14.7 (psia) 
6) wt = .22 
i) Ne = 02 5 


8) noscrw = l 
2) Dag eo Oe (re) 
10) Z= 6 
11) promat = 5 (stainless steel, see Table (II) ) 
meee D = 22.0 (ft) 
Two problems are examined. Problem 1 specifies the following 
additional parameters: 


me 6td 22 


By) =P 21292.6 (hp) 


E 
Some ¥ = 24.24 (knots) 


mee 00606. 75 (ft—ibt) 


eS 
© 
WY) 
| 


e) «UN 


OSs) (5pm) 
Problem 2 specifies the same parameters as: 
meeeca = .1l/l 
2) Po = 7630.0 (hp) 
See VY = 23.0 (knots) 
4) o> ra060G"75 (fLt=lbe ) 
5) N» = 105 (xm) 
All of the above are initialized in the input section 
(ICALC = 1) of similar versions of SUBROUTINE ANALIZ. There- 
fore, only one version is included in Appendix I. 


ae constraints for Gy (X3) and Gy 5 (X3) are evaluated by 


SUBROUTINE BLPOW3 which appears in the execution section of 
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SUBROUTINE ANALIZ. Also, note that SUBROUTINE DICNUA has 
been deleted from the execution section, while SUBROUTINE 
WGTCAL has been added. 


lr, Programming Details 


All twelve constraints are evaluated (NCON = 12). 


lower 
1 


upper 


The upper (X3, ) and lower (X3 


) limits on the design 


, * ° 
variables P/D, AL/Ap and (t /C) JER are set to be: 


4 = P/D < 1.4 
4 = AD/A, Se 
oO <4 t-7 CG) 75R < Poe 


ese Upper and lower limits are specified in the COPES con- 
trol card deck on card image F under respective fields VUB 
and VLB. The initial value for each design variable (X3.) 
is also assigned on card image F under the field labeled xX. 
menwrist of card images in Appendix J lists all of the COPES 
control cards used for both problems. These cards also 
specify the locations of the design variables in the common 
block GLOBCM (see Table (III)) as well as the locations of 
the constraints and their boundaries. Further details on 
the COPES control card requirements and the format of each 
Card are contained in Reference [7]. 
2. Results 

The outputs from the optimization/analysis, performed 

by COPES/CONMIN, are listed in Appendix K. Results of both 


problems are tabulated in Table (VI). 
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mee DISCUSSION 

Table (VI) presents the results of problems 1 and 2 
along with relevant information from Vassilopoulos' Deen 
Problem 1 attempted to "match" a Wageningen propeller at the 
design point found by Vassilopoulos. The first COPES/CONMIN 
printout in Appendix K indicates that the "match" was achieved 


Semeeeeequal to 21.168.1 (hp) and P, equal to 28,150.0 (hp) 


E D 
HOr, Qo = 1500607 (ft-lbf) and Np = 105 (rpm)). These values 
are judged to be close enough to the "Given" values in Table 
m1). 

It is apparent that the Wageningen propeller does not 
require all of the 30,000 (hp) of delivered horsepower. The 


propeller characteristics (i.e., J, K K.~ and No) for problem 


oc 
1 compare very well to Vassilopoulos' values. The expanded 
area ratios (Ap/Ad) are, also, very similar. Of course, the 
obvious difference is the blade weight (bldwt). The Wageningen 
propeller blade is over five thousand pounds heavier. Does 
this make sense for a minimum blade weight? 

The answer is yes. 
All one has to do is consider the values of (t*/c) for problem 
1 and Vassilopoulos' design. Vassilopoulos' "constant stress" 
blade was designed to "absorb" stress up to the allowable 
Seemegm Limit of 5,400 (psi) (for stainless steel) all along 
the entire propeller radius (R). Table (12) in Reference 
[18] gives further details. The Wageningen propeller blade, 


however, represents an "older" type of blade which was de- 


Signed with a linear blade section maximum thickness (t*) 
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distribution. Consequently, it was "overdesigned" for strength 
Beyond the 3/10--4/10 radius (1.e., .3R--.4R) and contains 
excess material. A heavier blade, therefore, results. Note, 
also, that the optimizer did not drive the value of (t*/C) sep 
to the minimum acceptable value, (t*/C) sep min’ 

The results of problem 2 show the effect on blade weight 
(bldwt) for a Wageningen propeller when the hull's powering 
requirements (1.e., Po at V) have been reduced. The weight 
reduction of 2000 pounds is significant. The complete re- 
Suets are listed in the second COPES/CONMIN printout in 


Pependix K. 
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GROUP 


Given 


Design 
Variable 
Specified 


Design 
Variables 


Minimize 


Maximize 


Restric- 
tions 


Other 


Design Case No. 


ur Sa 


P/D 
AWAD 


BSc) =o5R 


bldwt 


= 
O io 


° 


* 
P/D varies with R 


PABiE Vi 


3--Results 
VASSILOPOULOS PROBLEM 
ile 2 

21292.6 21292-6 7630.0 
24.24 24.24 23720 
PS00607 1500607 1500607 

105 EOS 10S 
22.0 22,56 22.0 
* i, Leis 1.0906 
. 7/67 . 7944 tae 
.040 .0794 .0681L 
— 12842.6 10464.7 

691 — _— 
__ foDU5 B22 
= .0691 .068l 
woo. . 8290 . 1866 
aoa? ~2349 .2063 
.0478 .0448 Peo 7Z 
-— .6915 . 6950 
Wo? «2 — — 
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X. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The general purpose non-linear optimizer/synthesizer 
COPES/CONMIN has been successfully applied to three typical 
preliminary ship design propeller selection problems in which 
the Wageningen B-Screw Series is used. The formulation and 
programming of each required analysis code (i.e., SUBROUTINE 
ANALIZ) have been made as general as possible to allow the 
designer a broad variety of solution options for solving 
propeller selection problems which can be classified under 
any of the three Design Cases that were considered. The 
analysis codes have been "modularized" to the extent that 
methodical series data from other propeller series, which 
are available in the polynomial expression format of the 
B-Screw Series, can be easily adapted for powering analysis 
utilizing deSign optimization methods. 

Further flexibility in the solution to the propeller 
selection problem has been achieved by using COPES/CONMIN 
as the optimizer/synthesizer. The designer has now been 
afforded the additional capability of specifying the design 
variables, the objective functions and the constraints of his 
choice. By solving propeller selection problems in the way 
presented in this thesis, repetitive problem formulation and 
coding have been eliminated. 

There are other advantages to solving propeller selection 


problems specifically with COPES/CONMIN which have not been 
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directly addressed in this study. As stated in Chapter II, 
COPES/CONMIN is capable of performing optimization analyses, 
sensitivity studies, optimum sensitivity studies and optimi- 
zation using approximation technigues. The designer, there- 
fore, can select and perform any of these options, using the 
Seemanalysis codes which have been presented in this thesis. 
While the utilization of a general purpose non-linear 
optimizer in solving propeller selection problems allows the 
designer greater flexibility in the selection procedure, there 
1s one important limitation that should be stressed at this 
point. This concerns the question whether or not the solution 
vector, determined by the optimizer, is a "global" optimum. 
As stated in Chapter II, COPES/CONMIN assures that, if a 
feasible solution vector is found, it is, at least, a "local" 
Minimum (or maximum). This implies that, for two different 
initial design vectors which are specified in the COPES Con- 
trol Card deck on card image F, the same optimum solution 
May not be determined by the optimizer. Both solutions would 
correspond to minimums (or maximums) of the objective function 
and are, therefore, correct. But, does one or the other 
correspond to the minimum (or maximum) of the entire vector 
design space, i.e., the "global" optimum? For the moment, 
at least, there is no definitive answer to the question. 
Despite this uncertainty, progress in the field of design 
Optimization continues to be made. Current developments 


[Ref. 47] will soon allow the designer to have a choice in 
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selecting a specific optimization algorithm froma "library" 
of proven optimization programs which employ the latest 
state-of-theart numerical techniques. Again, uSing one 


analysis code, the designer will be able to generate any 


number of optimized solutions for the problem under study. 


B. RECOMMENDATIONS 

For future consideration, it 1S recommended that the 
automated design and trade-off capability, provided by a 
general purpose non-linear optimizer/synthesizer such as 
COPES/CONMIN, be applied to the more difficult problem of 
propeller design. 

As pointed out in Chapter I, the use of the Wageningen 
B-Screw Series represents a "selection" procedure rather 
than a "design" process. Today, analytical propeller design 
m@eeceaures, utilizing lifting line and lifting surface 
theory, are becoming increasingly popular among propeller 
designers. The propeller design, which results from the 
utilization of these analytical methods, is, unquestionably, 
more efficient than the standard series propeller. However, 
these methods require consideration of many more deSign varia- 
bles in the design process. This appears to be a natural 
application for the use of a general purpose non-linear 
optimizer/synthesizer. 

Here, an analysis code, much larger than those which 
have been presented in this study, could be developed which 


would incorporate the lifting line/lifting surface theory 
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for the determination of the propeller performance character- 
istics, the local cavitation numbers and also the calculation 
of the pressure distributions over the blade. These pres~ 
Sure distributions would be utilized in the strength analysis 
of the blade. This analysis would utilize the finite element 
technique on an appropriately generated mesh model of the 
blade. Having defined the steps for this design procedure 

in the analysis code, the propeller designer now "couples" 
his analysis to the optimizer/synthesizer for determination 
of the optimum design. A massive amount of computer storage 


would certainly be required, but this concept is feasible 


and, in the author's view, is worthy of future consideration. 


C. A FINAL NOTE 

Poe conclusion, this thesis has demonstrated, in effect, 
another interesting application of the method of design opti- 
mization. The author, in no way, wishes to leave the reader 
with the impression that the techniques of design optimi- 
zation are the "be all-~end all" for engineering analysis. 
DeSign optimization techniques are useful and powerful tools 
that stand to relieve the engineer of the mundane tasks of 
numerical calculations and subsequent graphic plotting. 
But, they are just tools. In the final "analysis", good 


engineering judgment is paramount in their application and 


use. 
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APPENDIX A 


FORTRAN VARIABLE CROSS REFERENCE LIST 


Symbol 
Ap/Ag 


(AB/Ad) min 
blidwt 


 75R 


D 


ra 


Dog 


Pwatvap 


Patm 


promat 
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Fortran Variable 


AEDVAO 
AEAOMN 
WEIGHT 
CyoR 
DIA 
DIALIM 
reg 

RJ 

KQ 

Ke 
NOSCRW 
N 

EE 

Es 
EDLY) 
PWATVA 
PATM 
PROMAT 
QS 
Reo 
SC 


TD 





Symbol 
Temp 


(t*/C) FSR 
V (ft/sec) 
V (knots) 


wt 


APPENDIX A (CONT.) 
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Fortran Variable 
TEMP 
TEVOR 
V 
VK 


WT 


ETAO 
ETARR 
WATNU 


WATRO 
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OEFINITION 
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SUBROUTINE: 
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~<I 2 


= J 


UPELLER 


A. 
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